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The multivariate Tutte polynomial (alias Potts model) 
for graphs and matroids 

Alan D. Sokal 

Abstract 

The multivariate Tutte polynomial (known to physicists as the Potts-model 
partition function) can be defined on an arbitrary finite graph G , or more gen¬ 
erally on an arbitrary matroid M, and encodes much important combinatorial 
information about the graph (indeed, in the matroid case it encodes the full 
structure of the matroid). It contains as a special case the familiar two-variable 
Tutte polynomial — and therefore also its one-variable specializations such as 
the chromatic polynomial, the flow polynomial and the reliability polynomial — 
but is considerably more flexible. I begin by giving an introduction to all these 
problems, stressing the advantages of working with the multivariate version. 

I then discuss some questions concerning the complex zeros of the multivariate 
Tutte polynomial, along with their physical interpretations in statistical me¬ 
chanics (in connection with the Yang-Lee approach to phase transitions) and 
electrical circuit theory. Along the way I mention numerous open problems. 

This survey is intended to be understandable to mathematicians with no prior 
knowledge of physics. 

1 Introduction 

Let G = (V, E ) be a finite undirected graph with vertex set V and edge set E. 1 
The multivariate Tutte polynomial of G is, by definition, the polynomial 

S G (9,V) = (1-1) 

ACE ecA 

where q and v = {f e }ees are commuting indeterminates, and k(A) denotes the 
number of connected components in the subgraph (I/, A). [It is sometimes convenient 
to consider instead 

Z G (q,v) = q-Wz G (q,v) = ^ q HA)-\V\ JJ Ve ; (L2) 

ACE e€ A 

which is a polynomial in q^ 1 and {u e }.] From a combinatorial point of view, Zq is 
simply the multivariate generating polynomial that enumerates the spanning sub¬ 
graphs of G according to their precise edge content (with weight v e for the edge e) 
and their number of connected components (with weight q for each component). 
As we shall see, Zq encodes a vast amount of combinatorial information about the 
graph G , and contains many other well-known graph polynomials as special cases. 
I shall most often take an analytic point of view, and treat q and {v e } as real or 
complex variables. 


1 In this paper a 
otherwise. 


“graph” is allowed to have loops and/or multiple edges unless explicitly stated 
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In statistical physics, Z G is known as the partition function of the q-state Potts 
model [102,140,141], for reasons to be explained in the next section. 2 The Potts 
model — along with its close relative, the Fortuin-Kasteleyn random-cluster model 
[58,62,76] - plays an important role in the theory of phase transitions and critical 
phenomena [5,44,96], and physicists have developed a significant amount of infor¬ 
mation (both rigorous and non-rigorous) about its properties. As will be explained 
in Section 5, one way of analyzing phase transitions is to study the zeros of Zc(q, v) 
when q and/or {u e } are treated as complex variables. 

If we set all the edge weights v e equal to the same value v, we obtain a two- 
variable polynomial Zc{q,v) that is essentially equivalent to the standard Tutte 
polynomial Tq(x, y ) [see Section 2.5 for details]. But the main message of this paper 
is that it is often useful to consider the multivariate polynomial Z G (q,v), even if 
one is ultimately interested in a particular two-variable or one-variable specialization. 
For instance, Zo(q,v) is multiaffine in the variables v (i.e., of degree 1 in each v e 
separately); often a multiaffine polynomial in many variables is easier to handle 
than a general polynomial in a single variable (e.g., it may permit simple proofs by 
induction on the number of variables). Furthermore, many natural operations on 
graphs, such as the reduction of edges in series or parallel, lead out of the class of 
“all v e equal”. I shall illustrate the advantages of this “multivariate ideology” by 
showing several instances in which the multivariate extension of a single-variable 
result is not only vastly more powerful but also much easier to prove: indeed, in 
one case, a 20-page proof is reduced to a few lines. 

All of these considerations can be extended from graphs to matroids 3 : if M is a 
matroid with ground set E and rank function tm , its multivariate Tutte polynomial 
is defined by 

Z M (q,v) = Y, IK, (1-3) 

ACE eCA 

which is a polynomial in q~ l and {v e }. This extends the graph definition in the 
sense that if G is a graph and M(G ) is its cycle matroid, then 

Z M (G) (<?, v ) = Z G (q,v) (1.4) 

[because r M ^(A) = |V| — k(A)]. Since a matroid is completely determined by its 
rank function, Zm is simply an algebraic encoding of all the information about the 
matroid M. Moreover, my earlier statement that Zq encodes “a vast amount” of 
information about the graph G can now be made more precise: Zq encodes the 
number of vertices V together with all the information about G that is contained 
in its cycle matroid M(G ) [and no other information], 

I am now convinced that matroids are the natural category for studying the 
multivariate Tutte polynomial. Of course, results that hold for graphic matroids 

2 The Potts model [102] was invented in the early 1950s by Potts’ thesis advisor Domb: see [47]. 
The q — 2 case, known as the Ising model [67], was invented in 1920 by Ising’s thesis advisor 
Lenz [84]: see [20, 78,127] for a fascinating history. (I hasten to add that these are the only two 
cases I know of where the thesis advisor’s invention was named after the graduate student, rather 
than the other way around.) The q = 4 case, which is a special case of the Ashkin-Teller model, 
was invented in 1943 by Ashkin and Teller [3]. For a review of the Potts model, see [140,141]. 

3 See Oxley [98] for an excellent introduction to matroid theory. 
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may or may not hold for larger classes of matroids. 


2 The multivariate Tutte polynomial for graphs, and its specializations 

First, a trivial but useful observation: Starting from the definition 

Z G (q,v) = ^ H v * ( 2 - 1 ) 

ACE eCA 

and using the simple relation 


k(A) = \V\ - \A\ + c(A) (2.2) 

where c(A) is the cyclomatic number (i.e., number of linearly independent cycles) 
of the graph (V. A), we can rewrite Zq in the alternate form 

z G (q,v) = q\ y \ E^ (A) n-- ( 2 - 3 ) 

ACE eCA q 

Example 2.1 Trees. For any tree T = (V.E) we have 

Z T {q,v) = gfj(g + v e ). (2.4) 

e£E 

This follows immediately from (2.3), using the fact that c(A) = 0 for all A. □ 


Example 2.2 Cycles. Let G = (V,E) be the cycle C n . Then 

Zc n {q,v) = JJ(g + v e ) + (q - 1) JJ v e . (2.5) 

eCE eCE 

This follows immediately from (2.3), using the fact that c(A) = 0 for all A C E and 
c(E) = 1. □ 

Let us now consider some of the polynomials that can be obtained from Z G (q, v) 
by specializing the value of q: 


2.1 q = 1 

When q = 1, the multivariate Tutte polynomial becomes trivial: 

Zg{ !,v) = J|(l+u e ). (2.6) 

eCE 


2.2 q = 1, 2, 3,... (g-state Potts model) and the chromatic polynomial 

Let q be a positive integer; then the q-state Potts-model partition function for 
the graph G is defined by 


Z, 


Potts 

G 


(9,v) = 


e n 


1 T V e S(a Xl ( e ), &x 2 (e)) 


a: V->{1,2eCE L 


(2.7) 
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Here the sum runs over all maps a: V —> {1, 2 ,,q}, and we sometimes write a x 
as a synonym for cr(x); the 5 is the Kronecker delta 


6(a, b) 


1 if a = b 
0 if a ^ b 


( 2 . 8 ) 


and xi(e),X 2 (e) € V are the two endpoints of the edge e (in arbitrary order). We 
usually consider v e to be a real or complex variable. 

In statistical physics, this formula arises as follows: In the Potts model [102,140, 
141], an “atom” (or “spin”) at the site x £ V can exist in any one of q different 
states. The energy of a configuration is the sum, over all edges e € E, of 0 if the 
spins at the two endpoints of that edge are unequal and —J e if they are equal. 
The Boltzmann weight of a configuration is then e~^ H , where H is the energy of 
the configuration and /3 > 0 is the inverse temperature. The partition function is 
the sum, over all configurations, of their Boltzmann weights. Clearly this is just a 
rephrasing of (2.7), with v e = e^ Je — 1. A coupling J e (or v e ) is called ferromagnetic 
if Je > 0 (v e > 0), as it is then favored for adjacent spins to take the same value; 
antiferromagnetic if —oo < J e < 0 (—1 < v e < 0), as it is then favored for adjacent 
spins to take different values; and unphysical if v e £ [—1, oo), as the weights are then 
no longer nonnegative. 

It is far from obvious that Z^T tts (q, v), which is defined separately for each 
positive integer q, is in fact the restriction to q £ Z + of a polynomial in q. But this 
is in fact the case, and indeed we have: 


Theorem 2.3 (Fortuin—Kasteleyn representation of the Potts model) 

For integer q > 1, 

Z^ tts (q,v) = Z G (q,v) . (2.9) 

That is, the Potts-model partition function is simply the specialization of the multi¬ 
variate Tutte polynomial to q £ Z+. 

Proof In (2.7), expand out the product over e £ E, and let A C E be the set 
of edges for which the term v e S(a Xl r e \, cr X2 ( e )) is taken. Now perform the sum over 
configurations {<J x } xe v'■ hr each component of the subgraph (V, A) the color a x must 
be constant, and there are no other constraints. Therefore, 

Z G (q,v) = V HA) IK’ ( 2 T 0 ) 

ACE eCA 

as was to be proved. □ 


The subgraph expansion (2.10) of the Potts-model partition function was dis¬ 
covered in the late 1960s by Fortuin and Kasteleyn [58, 76]. Please note that if 
q > 0 and v e > 0 for all e, then the weights in (2.10) are nonnegative, and so can be 
interpreted probabilistically (after normalization by Z G ). The resulting probability 
measure on 2 E is called the Fortuin-Kasteleyn random-cluster model [58,62,76] for 

4 In the special case v e = — 1 it was discovered many decades earlier by Birklioff [8] and Whitney 
[136]; see also Tutte [130,131]. 
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the graph G. Furthermore, for q G Z+ a natural joint probability measure on spin 
states and edge occupations [i.e., a coupling of (2.7) and (2.10)] has been given by 
Edwards and Sokal [53]; it underlies the Swendsen-Wang [128] Monte Carlo algo¬ 
rithm for the ferromagnetic Potts model, and is also useful in the rigorous analysis 
of the Potts model. 

Another important special case arises when v e = — 1 for all edges e: then Z^ otts 
gives weight 1 to each proper coloring and weight 0 to each improper coloring, and 
so counts the proper colorings. It follows from Theorem 2.3 that the number of 
proper (/-colorings of G is in fact the restriction to q G of a polynomial in q, 
called the chromatic polynomial Pg(q) = Zc(q, — 1). The chromatic polynomial 
thus corresponds to the zero-temperature (/? —> +oo) limit of the antiferromagnetic 
(J e < 0) Potts model. Many properties of the chromatic polynomial in fact extend 
to the entire antiferromagnetic region (i.e., — 1 < v e < 0 for all e). 

2.3 q —» 0 limits 

Let us now consider the different ways in which a meaningful q —> 0 limit can be 
taken in the (/-state Potts model. 

The simplest limit is to take q —> 0 with fixed v. From the definition (2.1) we 
see that this selects out the subgraphs A C E having the smallest possible number 
of connected components; the minimum achievable value is of course k(G ) itself (= 
1 in case G is connected). We therefore have 

lim [q ~ k ^Z G (q,v) = C G (v), (2.11) 

q^O 

where 

CgM = J2 U y e ( 2 - 12 ) 

ACE eeA 
k(A) = k(G) 

is the generating polynomial of “maximally connected spanning subgraphs” (= con¬ 
nected spanning subgraphs in case G is connected). 

A different limit can be obtained by taking q —► 0 with fixed values of w = 
v/q. From (2.3) we see that this selects out the subgraphs A C E having the 
smallest possible cyclomatic number; the minimum achievable value is of course 0. 
We therefore have [126,139] 

limqH 1 \Z G {q, <?w) = F G ( w) 

q^O 

where 

f g( w) = II We 

ACE eeA 
c(A) = 0 

is the generating polynomial of spanning forests. 

Finally, suppose that in C G (v) we replace each edge weight v e by Xv e and then 
take A —► 0. This obviously selects out, from among the maximally connected 
spanning subgraphs, those having the fewest edges: these are precisely the maximal 


(2.13) 

(2.14) 
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spanning forests (= spanning trees in case G is connected), and they all have exactly 
\V\ — k(G) edges. Hence 


lim A fc(G )-l y IC' G (Av) = T G (v) , (2.15) 

A—>0 

where 

T G (y) = Y II* (2-16) 

ACE eeA 
k(A ) = k(G) 
c(A) = 0 

is the generating polynomial of maximal spanning forests. 5 Alternatively, suppose 
that in Fq( w) we replace each edge weight w e by \w e and then take A —> oo. This 
obviously selects out, from among the spanning forests, those having the greatest 
number of edges: these are once again the maximal spanning forests. Hence 

lim A fc ( G )“^.F G (Aw) = T G (w) . (2-17) 

A—xx) 

In summary, we have the following scheme for the q —> 0 limits of the multivariate 
Tutte polynomial: 


C'g(v) 


q —> 0, v fixed 


Zg(<-1;V) 


q —> 0, w = v/q fixed 


Fg( w) 


r infinitesimal 


T G (v orw) (2-18) 


r infinite 


Finally, maximal spanning forests (= spanning trees in case G is connected) can 
also be obtained directly from Z G (q,v) by a one-step process in which the limit 
q 0 is taken at fixed x = v/q a , where 0 < a < 1 [58,63,126,139]. Indeed, simple 
manipulation of (2.1) and (2.2) yields 

Z G {q,q a x) = q a \ y \ Y q ac ^ + ^- a ^ A ^ ]]_ x e . (2.19) 

ACE eCA 

The quantity ac(A) + (1 — a)k(A) is minimized on (and only on) maximal spanning 
forests, where it takes the value (1 — a)k(G). Hence 

limq-“ |v| - (1 -" )fc(G) Z G (g,g Q x) = T G (x) . (2.20) 

q-^ o 


5 I trust that there will be no confusion between the generating polynomial Tg(v) and the Tutte 
polynomial Tc(x,y). I have used here the letter T because in the most important applications the 
graph G is connected, so that Tg(v) is the generating polynomial of spanning trees. 
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2.4 The multivariate flow polynomial 

Let G = (V, E) be a finite undirected graph, and let T be a finite abelian group 
of order q = |T|. Let us choose arbitrarily an orientation for each edge of G (all 
subsequent results will be independent of this choice). A T-flow on G is a mapping 
ip: E —> r that satisfies current conservation at every vertex. A T-flow ip on G is 
said to be nowhere-zero in case ip(e) 0 for all e G E. 

Let Fq(T) be the number of nowhere-zero T-flows on G. 6 It is a well-known 
(but at first sight quite surprising) fact that Tb(T) depends only on the order of 
T and not on any other aspect of the structure of T; so we can write it as Fc{q). 
Furthermore, it turns out that F G (q) is the restriction to Z + of a polynomial in q, 
called the flow polynomial of G and also written Fc(q). 

In order to prove these facts (and much more), let us generalize F(-;(T) as follows: 
Assign to each edge e € E a weight u e , and write u = { u e } e £E■ Then define 


f g { r,u) 


E 

r-flows i/j 


n 


1 + u e 6(ip(e), 0 ) , 


e£E 


( 2 . 21 ) 


where the sum runs over all T-flows ip on G. If we take u e = — 1 for all e, this 
reduces to our friend Fq(T). We now assert that Fg(T,u) depends only on q = |T| 
- so we can write it as F G (q, u) — and that Fc(q, u) is in fact the restriction to 
q G Z + of a polynomial in q and {u e }, which we shall call the multivariate flow 
polynomial. Indeed, the multivariate flow polynomial is nothing other than the 
multivariate Tutte polynomial in disguise: 


Theorem 2.4 Let G = (V, E) be a finite undirected graph, and let T be a finite 
abelian group of order q = |T|. Then 


F G (F,u) 



( 2 . 22 ) 


(where division is of course understood edge-wise) 


and in particular 


F G (T ) = q~ l y ' (-I)'*' Z G (q, -q ) . (2.23) 


Consequently, Tb(T) and Fq(T, u) depend only on q = |T| and not on any other 
aspect of the structure ofT. 


Proof Note first that the total number of T-flows on a graph G is |T| C ( G ), where 
c{G) is the cyclomatic number of G. Moreover, the T-flows that vanish on a subset 
A C E are in one-to-one correspondence with the T-flows on the subgraph (V, E\A), 


6 I trust that there will be no confusion between the flow polynomial Tg(T) [and its relatives to 
be defined below] and the spanning-forest polynomial Fg( w). 
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so their number is |r| c ^^\ Expanding the product in (2.21), we find 


F G (X, u) 


by comparison with (2.3). 



c(E\A) JJ 

U e 

(2.24a) 

ACE 

eCA 



(n 


q< A "> TT 1 

(2.24b) 

\eeE 

/ A'CE 

e£A' e 


q- m 

fn«.) 

Zc(q,q/ u) 

(2.24c) 


\e£E / 






□ 


Theorem 2.4 is in fact a special case of a more general identity [31,49,113,142] 
relating “spin models” taking values in an abelian group T to “flow models” taking 
values in the dual group T*: here the Boltzmann weight W e (a x — a y ) for an edge 
e = xy in the spin model maps onto a weight ll / e (V’(e)) in the flow model, where 
" denotes Fourier transformation. In the special case of the Potts model, we have 
W e = 1 + v e 5 , so that W e = v e l + q5. Setting u e = q/v e and relabelling T* as T, we 
obtain (2.22). 


2.5 Comparison to the standard Tutte polynomial 

The Tutte polynomial Tc(x,y ) is conventionally defined by [135, p. 45] [28, 
pp. 124-127] 


{X - 1 )r(^)-r(A) ( y _ 1 j|A|-r(A) 

ACE 

(2.25a) 

^2(x- l) fc ( A )- fc (S) (y _ 1 )|A|+fc(A)-|V| 

ACE 

(2.25b) 


where r(A ) = \V\ — k(A) is the rank of the set A in the cycle matroid M(G). [Note 
also that |A| — r{A) is the cyclomatic number c(A).] Comparison with (1.1) yields 


Tc(x,y ) = {x 1) HE) {y- 1) |v 1 Z G ((x - l)(y - 1), y - l) . 


(2.26) 


In other words, Tc(x,y) and Za(q,v) are essentially equivalent under the change of 
variables 


x = 

1 + q/v 

(2.27a) 

y = 

1+v 

(2.27b) 

q = 

(■x - 1 ){y - 1) 

(2.27c) 

V = 

y- 1 

(2.27d) 


One advantage of the Tutte notation is that it allows a slightly smoother treat¬ 
ment of the q —> 0 limit, in which the limiting processes employed in Section 2.3 
can be replaced by simple evaluation at x = 1 or y = 1. Thus, the univariate 
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maximally-connected-spanning-subgraph and spanning-forest polynomials are given, 
respectively, by 


C G {v) = r G (l, 1 + v) (2.28) 

F g (w) = w lvl ~ k{E) T G ( l + l/w,l) (2.29) 

A second advantage of the Tutte notation is that duality [cf. (4.7) below] takes the 
simple form 

T G * (x, y) = T G (y, x) . (2.30) 

But the Tutte notation also has a severe disadvantage: the use of the variables x 
and y conceals the fact that the particular combinations q and v play very different 
roles', q is a global variable, while v can be given separate values v e on each edge. 
In particular, the Tutte notation makes it impossible to assign unequal weights v e 
to the edges. I therefore recommend use of the notation Z G (q,v) whenever the 
“multivariate ideology” is potentially of use. 

Let me conclude by observing that numerous specific evaluations of the Tutte 
polynomial have been given combinatorial interpretations, as counting some set of 
objects associated to the graph G (see e.g. [28,135]). It would be an interesting 
project to seek to extend these counting problems to “counting with weights”, i.e., 
to obtain suitably defined univariate or multivariate generating polynomials for the 
objects in question as specializations of Z G (q,v) or Z G (q,v), respectively. A few 
examples of this are given in Sections 2.3 and 2.4 above. 

3 The multivariate Tutte polynomial for matroids 

Many (but not all) of the properties of the multivariate Tutte polynomial for 
graphs can be carried over to the matroid version (1.3). Where the treatment is 
essentially identical, I shall be brief. 

3.1 q = l 

When q = 1, the multivariate Tutte polynomial becomes trivial: 

Z M { l,v) = JJ(l + u e ). (3.1) 

e£E 


3.2 q = 1, 2, 3,... (generalized Potts models) 

For graphs, we saw in Section 2.2 that the multivariate Tutte polynomial Z G (q, v) 
for integer q > 1 counts (/-colorings of the vertices of G, with weights that depend on 
whether each edge e is properly or improperly colored [cf. (2.7)/(2.9)]. In particular, 
the chromatic polynomial P G (q ) = Z G (q, —1) counts proper (/-colorings of G. One 
consequence of this is that the nonnegative integer zeros of P G {q) are consecutive 
integers, namely 0,1, 2,... , y(G') — 1, where x(G) is the chromatic number of G 7 

'In fact, these are all the integer zeros of Pdq), as it is easy to show that Pg(q) has no negative 
real zeros (when G is loopless): see e.g. [106]. 
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Simple examples show that the latter property does not hold in general for 
the chromatic polynomial (also known as characteristic polynomial) of a nratroid, 
Pm(q) = — l). 8 For instance, for the uniform matroid U^.n we have 

n n 

Zu 2 , n {q,v) = g" 2 JJ(l + Vi) + (1 ~q~ 2 ) + {q- 1 ~ q^^Vi (3.2) 

i =1 i=l 

and hence 

Pu 2 , n {q) = q~ 2 (q- 1)(? + 1 -«), (3-3) 

whose roots are q — 1 and q = n — 1; so the consecutive-root property is violated 
for all n > 4. It therefore seems unlikely that Z^v/((/,v) can be given a “coloring” 
interpretation for arbitrary matroids M and arbitrary integers q > 1. Nevertheless, 
in some special cases such a representation can be given, as I would now like to 
sketch; this account is drawn from work in progress with Sergio Caracciolo and 
Andrea Sportiello [32]. 

The appropriate general context for this discussion is that of abelian-group-valued 
statistical-mechanics models , defined as follows [31]: Let V and E be finite index 
sets (for variables and interactions , respectively); let R be a ring with identity and 
let Ad be a unitary right A-module; let B = ( bi e )i & v,e£E be a matrix with elements 
in R ; let /i be Haar measure on Ad, considered as an additive abelian group; and for 
each e € E, let W e : Ad —> C be a function. We then define the partition function 

z r ,mm w) = f n^) ( 3 . 4 ) 

J eS E 'SeV' ' i&V 

where W = { W e } e ^E■ (If Ad is infinite, then the collection of functions W needs to 
satisfy suitable integrability conditions.) 

[We can alternatively let Ad be an arbitrary additive abelian group Q, and let R 
be some subring of the ring End(f/) of endomorphisms of Q. It is easy to see that 
the two formulations are equivalent.] 

Two important special cases are: 

1) R = Z and Ad is an arbitrary additive abelian group Q (on which Z acts in 
the obvious way). 

2) R is a field F, and Ad is a finite-dimensional vector space over F (which can 
be taken, without loss of generality, to be F k for some integer k > 0). 

A generalized Potts model is the special case of this setup in which Ad is finite, /j 
is normalized counting measure on Ad (i.e., each element of Ad gets weight l/|Ad|), 
and each W e is of the form 

W e (a) = 1 + v e 5(a, 0) (3.5) 

for some weights v = {v e } e& E- We then have 

zgswv) = \ m\-m y . n 

_ <7 : V—>JA e£E 

8 The chromatic (or characteristic) polynomial of a matroid M is ordinarily defined as Tm(<?) = 
q r( ' M )P M {q), so that it is a polynomial in nonnegative powers of q. This is a matter of taste. 


1 + V e S 


E 

iev 


Ojbjp, 0 


(3.6) 
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Note that if G = (V, E) is a graph, B is the directed vertex-edge incidence matrix 
for some orientation of G (hence R = Z), and A4 is a finite abelian group of order 
q = \M\, we have -^r°Jv( S b(v) = Z^fi tts (q, v ) = ^Z[fi tts (q, v). So the models (3.6) 
do indeed constitute a generalization of the Potts models (2.7). 

As far as I know, there does not exist, at this level of generality , a satisfactory 
analogue of the Fortuin-Kasteleyn representation (2.9). But at least in some special 
cases such a representation does exist [32]: 

Theorem 3.1 (Fortuin-Kasteleyn representation for generalized Potts models) 

(a) Let M he a regular matroid, let B be a totally unimodular matrix of integers 
representing M over Q, and let Q be a finite abelian group of order q = \Q\. 

Then 

^£°g tt s( v ) = Z M (q,v) ■ (3.7) 

(b) Let M be a matroid, let F be a finite field of order q = |F| [so that q is a 
prime power and F ~ GF(q)J, let B be a matrix representing M over F, and 
let A4 = F k for some integer k > 0. Then 

= Z M (q\v). (3.8) 

We recall that an integer matrix B is called totally unimodular if all its subdetermi¬ 
nants lie in {0,1, —1}. 

Proof In (3.6), expand out the product over e£f, and let A C E be the set of 
edges for which the term v e 5{■ • •) is taken. We then need to count the number of 
configurations cr = {cq}j e y satisfying Xaev a ihe = 0 for all e € A. In other words, 
we need to count the number of solutions cr € of the linear system aB = 0, 
where B is the submatrix of B consisting of the columns from A. 

Since both R = Z and R = F are principal ideal domains, the matrix B has a 
Smith normal form [2,19,89,95] 

S = PBQ = [^ ° ) (3.9) 

where P G GL(m,R), Q G GL(n,R) and D = diag(si,..., s r )\ here m = |P|, 
n = | A |, r is the determinantal rank of B, and si,..., s r 0 are the invariant factors 
of B. In case (a), since B is totally unimodular, we must have si,...,s r G {1, —1}; 
and by a suitable choice of P (or Q) we can arrange to have si = ... = s r = 1. 

In case (b), because R is a field, we can again arrange to have si = ... = s r = 1. 

Therefore, a vector cr G M v solves crB = 0 if and only if crP -1 is vanishing in its first 

r components. The number of such solutions is thus |A1|^ . On the other hand, 

since R is an integral domain, the determinantal rank r is also equal [2, pp. 205-209] 
to the column rank of B in the quotient field of R [namely, Q in case (a), and F 
itself in case (b)]; and this rank is by definition This proves (3.7)/(3.8). □ 

Theorem 3.1(a) implies, in particular, the known fact [28, Exercise 6.53(e,f)] that 
the chromatic polynomial of a regular matroid has the consecutive-root property. 

See also [40, Theorem III] [27, Theorem 12.4] for results related to Theorem 3.1(a), 
and [41, Theorem 16.1] [144, Theorem 7.6.1] for results related to Theorem 3.1(b). 
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3.3 q —► 0 limits 

The q —* 0 limits for the matroid multivariate Tutte polynomial (1.3) follow 
closely the pattern for graphs, with the obvious replacements 

maximally connected spanning subgraph —> spanning set 

spanning forest —> independent set 
maximal spanning forest —> basis 

Thus, the limit q —► 0 with fixed v selects out the subsets A C E of maximum 
rank, i.e. the spanning sets. We therefore have 

lim q rM ( E ' ) Z M {q,v) = S M {v) , (3.10) 

g^O 

where tm{E) is the rank of M, and 

-S'm(v) = ( 3 - n ) 

ACE ee A 
r M(A) = r M (E ) 

is the generating polynomial of spanning sets in M. 

Similarly, the limit q —> 0 with fixed w = v/q selects out the subsets ACE 
having the smallest value of \A\ — tm(A); the minimum achievable value is 0, and 
the sets attaining it are the independent sets. We therefore have 

lim Z M {q, </w) = Im( w) , (3.12) 

g^O 

where 

Im( w) = ^2 n- ( 3 - 13 ) 

ACE eeA 

4) = |^4| 

is the generating polynomial of independent sets in M. 

Finally, we have 

lim A~ ? ’ m ( £ ')»S'm(Av) = %(v) (3.14a) 

lim \~ rM ^ I M (Xw) = Bm( w) (3.14b) 

A—>oo 

where 

%(v) = ^ (3.15) 

ACE ee A 

|^ 4 | = r M (A) = rM{E) 

is the generating polynomial of bases in M. 
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In summary, we have the following scheme for the q —» 0 limits of the multivariate 
Tutte polynomial for matroids: 


Z M (q,v) 


Sm{v) 



q —> 0, w = v/q fixed 



Im{ w) 



v infinitesimal 


Bm(v or w) 



(3.16) 


3.4 A final remark 

Several alternative definitions of multivariate (or “weighted”) Tutte polynomials 
for graphs and/or matroids can be found in the literature. Most of these are es¬ 
sentially equivalent to Zm(< 7 , v) after suitable changes of variables. Exceptions can, 
however, be found in the work of Zaslavsky [145] and Bollobas-Riordan [9]; their re¬ 
sults are generalized and clarified in an illuminating recent paper of Ellis-Monaghan 
and Traldi [54]. 


4 Elementary identities 

I now wish to prove some elementary identities for the multivariate Tutte poly¬ 
nomial. There are two alternative approaches to proving such identities: one is to 
prove the identity directly for indeterminate (or complex) q, using the subgraph ex¬ 
pansion (1.1) or its generalization (1.3) to matroids; the other is to prove the identity 
first for positive integer q, using the spin/coloring representation (2.7), and then to 
extend it to general q by arguing that two polynomials (or rational functions) that 
coincide at infinitely many points must be equal. The latter approach is perhaps less 
elegant, but it is often simpler or more intuitive. However, only the former approach 
extends to matroids. 

One way to guess (albeit not to prove) an identity for matroids is to prove it 
first for graphs, and then translate it from Z G to Z G = q~ l v I Zq\ usually the latter 
identity carries over verbatim to matroids, mutatis mutandis. 

4.1 Disjoint unions and direct sums 

If G is the disjoint union of G\ and G* 2 , then trivially 

Z G (q,v) = Z Gl (q,v) Z G2 (q,v) . (4.1) 

That is, Zq “factorizes over components”. 

A slightly less trivial situation arises when G consists of subgraphs G\ and G -2 
joined at a single cut vertex x; in this case 

ZGi(q, v) Z G2 {q, v) 

q 


z G (q,v) 


(4.2) 
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This is easily seen from the subgraph expansion in the form (2.3). It is also easily seen 
from the coloring representation (2.7), by first fixing the color a x at the cut vertex 
and then summing over it; from this viewpoint, (4.2) reflects the S q permutation 
symmetry of the (/-state Potts model. 9 We summarize (4.2) by saying that Z G 
“factorizes over blocks” modulo a factor q. 

The identities (4.1) and (4.2) can be written in a unified form, by using Z G = 
I Zq\ in both cases we have 

Zg(q,v) = Z Gl (g, v) Z G2 (q, v) . (4.3) 

This, in turn, is a special case of the following obvious fact: if a matroid M is the 
direct sum of matroids and M 2 , then 

Z M {q,v) = Z Ml (q,v)Z M2 {q,v) . (4.4) 


4.2 Duality 

Suppose first that G = (V, E) is a connected planar graph. Consider any plane 
embedding of G, and let G* = ( [V*,E *) be the corresponding dual graph. There is 
a natural bijection between E and E* (namely, an edge e € E is identified with the 
unique edge e* G E* that it crosses), so we shall henceforth identify E* with E. Of 
course, the vertex set V* can be identified with the faces in the given embedding of 
G, so by Euler’s relation we have 

\V\ - \E\ + |C*| = 2 . (4.5) 


Consider now any subset ACE, and draw in G* the complementary set of edges 
(E\A). Simple topological arguments then yield the relations 

kc{A) = cg*{E\A ) + 1 (4.6a) 

k G *{E\A ) = cg(A) + 1 (4.6b) 


where k = components and c = cyclomatic number. [Note that (4.6a) and (4.6b) are 
equivalent, as a consequence of (4.5).] Substituting (4.6) into (2.1)—(2.3), we deduce 
the duality relation [140] 


Zc*{q, v) 



(4.7) 


In brief, duality takes v e 1 —> q/v e (and inserts some prefactors). Note that two 
applications of (4.7) lead us back to where we started, thanks to (4.5). In terms of 
Zq = q~ l v Zq , we have 


Zc*(q, v) 



(4.8a) 


(4.8b) 


9 More precisely, it reflects the symmetry of the spin model under a global transformation a y 1 —> 
ga y (simultaneously for all y £ V) that acts transitively on each single-spin space. 
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In the q —> 0 limit, we obtain the following special cases of (4.7): 


Cg * (v) 


F g * (w) 


7b* (v) 



(4.9) 


(4.10) 


(4.11) 


We can also relate the multivariate Tutte polynomial on G* to the multivariate flow 
polynomial on G: by (2.22) we have 


Z G *{q, v) = qF G (q, v) . 


(4.12) 


In particular, the chromatic polynomial on G* is essentially identical to the flow 
polynomial on G: 

P G *(q) = qP G {q) ■ (4.13) 

The chromatic polynomial and the flow polynomial are thus “dual” objects. 

Among graphs, only planar graphs have duals with good properties [43, sec¬ 
tion 4.6] [98, section 5.2]; one major advantage of considering matroids is that every 
matroid has a dual. The duality formula for the matroid multivariate Tutte poly¬ 
nomial is, not surprisingly, identical in form to (4.8): 

Z M *(q,v) = q ~ r M*( E ) | JJ Ve I Z M {q,q/v) (4.14a) 

\e€E ) 

= q rM(E) [W-Az M {q,qM ■ (4.14b) 

\e£E ^ / 

[Here vm{E) is the rank of M and tm*{E) is the rank of M*\ their sum is |i£|.] 
Indeed, (4.14) is an easy consequence of the definition (1.3) together with the formula 
for the rank function of a dual: 


r M *{A) = \A\ + r M {E\A ) - r M (E) . (4.15) 

4.3 Deletion-contraction identity 

If e G E, let G \ e denote the graph obtained from G by deleting the edge e, and 
let G/e denote the graph obtained from G \ e by contracting the two endpoints of e 
into a single vertex (please note that we retain in G/e any loops or multiple edges 
that may be formed as a result of the contraction). Then, for any e G E, we have 
the identity 

Z G (q,v) = Z G \ e (q,v^ e ) + v e Z G/e (q,v^ e ) (4.16) 

where v^ e = {vj}f & E\ e . This is easily seen either from the coloring representation 
(2.7) or the subgraph expansion (1.1). Please note that the deletion-contraction 
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identity (4.16) takes the same form regardless of whether e is a normal edge, a loop, 
or a bridge (in contrast to the situation for the usual Tutte polynomial Tq). Of 
course, if e is a loop, then G/e = G \ e, so we can also write Z G = (1 + v e )Z G \ e = 
(1 + v e )Z G / e . Similarly, if e is a bridge, then G \ e is the disjoint union of two 
subgraphs G\ and G *2 while G/e is obtained by joining G\ and G 2 at a cut vertex, 
so that Z G/e = Z G \ e /q and hence Z G = (1 + v e /q)Z G \ e = (q + n e )Z G/e . 

In terms of Z G = q~ ^ I Z G , the deletion-contraction identity takes the form 


Z G — Z G \ e + ~Z G / e 

if e is not a loop 

(4.17a) 

z g = Z G \ e + v e Z G / e 



O 

+ 
1—1 

II 



— (1 + v e )Z G / e 

if e is a loop 

(4.17b) 

as easily follows from (4.16) together with the counting of vertices in 
Not surprisingly, the deletion-contraction formula for matroids 
form to (4.17): 

G \ e and G/e. 
is identical in 

Sr 

iN 

£1 & 
+ 

NSJ 

II 

if e is not a loop 

(4.18a) 

Zm = -^M\e + Ve^M/e 



= (1 + v e )Z M \ e 



= (1 + v e )Z M / e 

if e is a loop 

(4.18b) 

This easily follows from the formulae for the rank function of a deletion or contrac¬ 
tion: if A C E \ e, then 

rM\e(A ) = r M {A ) 


(4.19a) 

(A\ fr M (A Ue)- 

r M /e(A) = < 

\r M (A U e) 

1 if e is not a loop 
if e is a loop 

(4.19b) 


Remark. Many treatments in the literature define the Tutte polynomial by the 
deletion-contraction identity (4.16) together with the initial condition Z G (q, v) = 
q ^ I for an edgeless graph (or the analogous thing for the standard Tutte polynomial 
T g ). But this approach has the disadvantage that one must prove that this Z G is 
well-defined, i.e. that the result does not depend on the order in which one applies 
(4.16) to the various edges. A much cleaner approach, it seems to me, is to define 
Z G (q,v) by the explicit formula (1.1), and then deduce the deletion-contraction 
identity as an immediate property. 

4.4 Parallel-reduction identity 

If G contains edges ei,e 2 connecting the same pair of vertices x, y, they can be 
replaced, without changing the value of Z, by a single edge e = xy with weight 


v, 


(1 + v ei ){\ + v e2 ) - 1 = v ei + v e2 + v ei v e2 . 


(4.20) 
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This is easily seen either from the coloring representation (2.7) or the subgraph 
expansion (1.1). More formally, we can identify the new edge e with (for instance) 
the old edge e\ after deletion of e 2 , and thus write 

Zg(q, v^ ei ,e 2 ,v ei ,v e2 ) = Z G \ e2 (q,v^ ei}e2 ,v ei +v e2 + v ei v e2 ) . (4.21) 

The parallel-reduction rule (v\,V 2 ) 1 —> u e fj with 1 + v e g = (1 + ui)(l + V 2 ) can be 
remembered by the mnemonic “1+u multiplies”. We write v\\\v 2 = (l+ui)(l+U 2 ) — 1. 

A virtually identical formula holds for matroids: if e\ and e 2 are parallel elements 
in a matroid M (i.e., form a two-elenrent circuit), then 

Z M (q,^e 1 ,e 2 , v e 1 ,Ve 2 ) = Z M \ e2 (q, V^ ei ,e 2 i V ei + V e2 + V ei V e2 ) . (4.22) 

The formula (4.22) also holds trivially if e\ and e 2 are both loops. 


4.5 Series-reduction identity 

We say that edges ei,e 2 £ E are in series (in the narrow sense) if there exist 
vertices x,y,z £ V with x j- y and y / z such that e\ connects x and y, e 2 connects 
y and z, and y has degree 2 in G. In this case the pair of edges ei, e 2 can be replaced, 
without changing the value of Z. by a single edge e = xz with weight 


v e 


V ei Ve 2 

q+v ei + v e2 


(4.23) 


provided that we then multiply Z by the prefactor q + v ei + v e2 . More formally, we 
can identify the new edge e with (for instance) the old edge e\ after contraction of 
e 2 , and thus write 


ZG(q,-Vjtei,e 2 ,V ei ,Ve 2 ) = (<? + V ei + V C2 ) Z G/e2 (q, v^ eij62 , V ei V e2 /(q + V ei + V C2 )) . 

(4.24) 

This identity can be derived from the coloring representation (2.7) by noting that 


<2 

E [i + ^ei 5(a x ,(j y )\ [l + v e2 5(a y ,a z )\ 

Oy = l 


= q + v ei + v e2 + v ei v e2 S(a x , a z ) 

^ei^e2 


= (q + v ei +v e2 ) 


1 + 


q + v ei + v e 


-6((t Xi a z 


(4.25a) 

(4.25b) 


Alternatively, it can be derived from the subgraph expansion (1.1) by considering 
the four possibilities for the edges ei and e 2 to be occupied or empty and analyzing 
the number of connected components thereby created. The series-reduction rule 
(r>i, V 2 ) >—> u e ff = V 1 V 2 /{q + v\ + V 2 ) can be remembered by the mnemonic “1 + q/v 
multiplies”: namely, 


1 + — = 

^efF 



i + ^ 


V2 


(4.26) 


We write v\ V 2 = v\V 2 /(q + v\ + U 2 ). 
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Consider now the more general situation in which {ei, e 2 } is a two-edge cut of G 
(not necessarily the cut associated with a degree-2 vertex y)\ we then say that e±, e 2 
are in series (in the wide sense). It turns out that the identity (4.24) still holds. To 
see this, let us prove the generalization of this identity to matroids. Let e± and e 2 
be series elements in a matroid M , i.e., suppose that {ei,e 2 } is a cocircuit. Then, 
for any A C E \ {e\, e 2 }, we have 

r M {A Uei) = r M (A U e 2 ) = r M (A) + l (4.27) 

(since the complement of a cocircuit is a hyperplane). A short calculation using 
(4.19b) with e = e 2 then yields 

Z M {q > v^ ei;62 , v ei , v e2 ) = -- — z M /e 2 {q^^ eit e 2 ,v ei v e J(q + v ei +v e2 )) . 

(4.28) 

The formula (4.28) also holds trivially if e\ and e 2 are both coloops. 

Please note that duality v t—> qjv interchanges the parallel-reduction rule (“1 + v 
multiplies”) with the series-reduction rule (“1 +q/v multiplies”). This is no accident, 
since we now see that parallel-reduction and series-reduction (in the wide sense) are 
indeed duals of each other: {ei,e 2 } is a circuit (resp. cocircuit) in M if and only if 
it is a cocircuit (resp. circuit) in M*. 

Remark. Using the series-reduction formula (4.24) together with the multi¬ 
variate approach, one can give a simple proof of (a generalization of) the Brown- 
Hickman [17] theorem on chromatic roots of large subdivisions [122, Appendix A], 
Likewise, using the parallel- and series-reduction formulae (4.21)/(4.24), one can 
give a simple proof of (a slight generalization of) Thomassen’s [129] construction of 
2-degenerate graphs with arbitrarily large real chromatic roots [122, Appendix B], 

4.6 Reduction formulae for 2-rooted subgraphs 

Let G = (V, E ) be a finite graph, and let x, y be distinct vertices of G. We define 
G/xy to be the graph in which x and y are contracted to a single vertex. (N.B.: 
If G contains one or more edges xy, then these edges are not deleted, but become 
loops in G/xy.) There is a canonical one-to-one correspondence between the edges 
of G and the edges of G/xy, for simplicity (though by slight abuse of notation) we 
denote an edge of G and the corresponding edge of G/xy by the same letter. In 
particular, we can apply a given set of edge weights {v e } e& E to both G and G/xy. 

Let us now define 

4r y) (q,v) = ]T q k(A > n«e (4.29) 

ACE eeA 

A connects x to y 

4^’(,,v) = £ /0> n „ e (4.30) 

ACE eeA 

A does not connect x to y 


= z, 


(x^y) 

G 


i,.v) + z<y*»>(,.v) 


From (1.1) we have trivially 
Zc(q,v) 


(4.31) 
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and almost as trivially 

Z G/xy(QiV) = Z^ y \q,v) + q~ 1 Z^ y \q,v) . (4.32) 

Let now q be an integer > 1, and define the restricted Potts-model partition 
function 


Z G°^(<2tv;cr 


= 


x •) w y 


E 

V\{x,y}-^ 


n 

{!....,(?} e£E 


1 "P V e S(^0’ Xl ( e ), &x2{e )) 


(4.33) 


where G x ,a y G {1,... , q}. We then have the following refinement of the Fortuin- 
Kasteleyn identity (2.9): 


Proposition 4.1 

Z G%%{^ v-,a x ,(Ty) = A GtXi y(q,v) + B GjXj y(q,v)8(a x ,a y ) (4.34) 

where 

A G , x , y (q,v) = q~ 2 Z^ y \q,v) (4.35a) 

B G , x , y (q,v) = q~ l Z^ y \q,\) (4.35b) 

are polynomials in q and {v e }, whose degrees in q are 

deg^4 = |P| — 2 (4.36a) 

degf? < \V\ — 1 — dist(x,y) (4.36b) 

where dist(a :,y) is the length of the shortest path from x to y using edges having 
v e 0 (if no such path exists, then B = 0). 


First Proof In (4.33), expand out the product over e G E, and let A C E be the 
set of edges for which the term v e 6(cr xi r e \, a X2 ( e )) is taken. Now perform the sum over 
configurations {&z} z ev\{x,y} : each connected component of the subgraph (V, A) 

the spin value a z must be constant. In particular, in each component containing 
x and/or y , the spins must all equal the specified value a x and/or a y ] in all other 
components, the spin value is free. Therefore, 


/7potts 
^G,x,y 


(q,v,(Tx, 


°y) 


q 2 Z^ y \q,v) if a x ^ a y 

<T 2 4^%, v ) + q- 1Z G~ v) (Q’ v ) if a * = a y 


(4.37) 


The claims about the degrees are easily verified. 


□ 


Second Proof For each positive integer q, the S q permutation symmetry of the 
Potts model implies that 

v;a x ,(Ty) = A + B8(a x ,a y ) (4.38) 

for some numbers A and B depending on G, x, y, q and v. But then, summing over 
a x , a y without and with the constraint a x = a y , we get 

= q 2 A + qB 
= qA + qB 


Z G 

Z G/xy 


(4.39a) 

(4.39b) 
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Hence 


A = 

B = 

by virtue of (4.31) and (4.32). 


ZG Z G/xy ^_ 2 ^(x^y) 

9(«-l) “ 9 ° 

(4.40a) 

qZc/xy ~ Zq ^_ 1 ^(x^y) 

9(9-1) - “ Z = 

(4.40b) 


□ 


Let us now consider inserting the 2-rooted graph (G, x, y) in place of an edge 
e* in some other graph H. From (4.34)/(4.35) we see that ( G,x,y ) then acts as a 
single edge with effective weight 


^eff — ^eS,G,x,y v) 


BG,x,y{q,v) 

A G,x,y(q,v) 


qZ£~ y) (q,v) 
Z^ y) (q,v) ’ 


(4.41) 


which is a rational function of q and { v e }; in addition, the partition function is 
multiplied by an overall prefactor A GtX ^ y (q,v). This follows from Proposition 4.1 
whenever q is an integer > 1; and the corresponding identity then holds for all q , 
because both sides are rational functions of q that agree at infinitely many points. 
It is also worth noting that the “transmissivity” t. e g = v e g/(q + u e ff) is given by the 
simple formula 

_ 

ZG(9.V) 



The most general version of this construction appears to be the following [138]: 
Let H = (V, E ) be a finite undirected graph, and let I? be a directed graph obtained 
by assigning an orientation to each edge of H. For each edge e € E, let G e = 
(V e , E e , x e , y e ) be a 2-rooted finite undirected graph (so that x e ,y e G V e with x e ^ y e ) 
equipped with edge weights {v^feeE e - We denote by G the family {G e } eG £, and we 
denote by H G the undirected graph obtained from H by replacing each edge e 6 1? 
with a copy of the corresponding graph G e , attaching x e to the tail of e and y e to 
the head. Its edge set is thus E = [+J E e (disjoint union). We then have: 

e£E 


Proposition 4.2 Let H = (V,E) and {G e } e ^E be as above. Suppose that 
Z Ge ,x e ,y e (q, { v e}e£E e i a x e , a i / e ) 


and define 


Then 




j e ,Xe,y e { c L { v e}e&E e ) b((T Xe , <Jy e ) , 

(4.43) 

{q, {Vehe Be) 

:( q,{vs}e€E e ) ' 

(4.44) 

e}e£E e )\ X Z H {q,{VeS,e}e£E) ■ 

(4.45) 


In particular, Zg G does not depend on the orientations of the edges of H. 
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4.7 Analogy with electrical circuit theory 

The identities discussed in Sections 4.4-4.6 are strongly reminiscent of analogous 
identities in electrical circuit theory: there are elementary formulae for the reduction 
of linear circuit elements placed in series or parallel; and more generally, any 2- 
terminal subnetwork consisting of linear passive circuit elements is equivalent (at a 
fixed frequency u) to some single “effective admittance”. But the relation of the 
multivariate Tutte polynomial to electrical circuit theory goes beyond mere analogy. 
For, as was discovered by Kirchhoff [77] in 1847 and will be reviewed in Section 6, 
linear electrical circuits are intimately related to spanning trees; and as was seen 
in (2.18), the spanning-tree polynomial 7 g(v) arises from the double limit of the 
multivariate Tutte polynomial in which q —> 0 and v is infinitesimal. 

Look back, for instance, at the Tutte-Potts parallel law v\\\v 2 = V\+V 2 +V\V 2 ', it is 
a nonlinear generalization of the familiar law v\ \\v 2 = V 1 + V 2 for combining electrical 
conductances in parallel, and reduces to it when v is infinitesimal. Likewise, the 
Tutte-Potts series law v\ \x\q V 2 = v\V 2 /{q + v\ + V 2 ) is a nonlinear generalization of 
the familiar law = ^ 1 ^ 2 /("^l + ^ 2 ) f° r combining electrical conductances in 

series, and reduces to it when q —> 0. 

It thus makes sense to ask which results of electrical circuit theory can be general¬ 
ized from the spanning-tree polynomial Tq(v) to the connected-spanning-subgraph 
polynomial C'g(v), the spanning-forest polynomial Fq{ w), or even the full multi¬ 
variate Tutte polynomial Zc(q,v). 

5 Complex zeros of Z G : Why should we care? 

Since Z G (q, v) is a polynomial in q and {v e }, it is mathematically natural to 
inquire about the zeros of Z G (q,v) as a function of q or {v e } or both. Here we can 
treat q and {u e } as either real or complex variables. 

The real zeros of the chromatic polynomial have been extensively studied, and 
to a lesser extent this is true also for the flow polynomial; see [70] for an excellent 
recent survey. Many of these results for the chromatic polynomial extend to the 
multivariate Tutte polynomial Z G (q, v) in the antiferromagnetic regime (i.e., — 1 < 
v e < 0 for all e) [74], 10 

The complex zeros of the partition function play a very important role in statis¬ 
tical mechanics, as was shown by Yang and Lee [143] in 1952. This arises as follows: 
Statistical physicists are interested in phase transitions, namely in points where one 
or more physical quantities (e.g. the energy or the magnetization) depend nonan- 
alytically (in many cases even discontinuously) on one or more control parameters 
(e.g. the temperature or the magnetic field). Now, such nonanalyticity is mani¬ 
festly impossible in (1.1)/(2.7) — or in any other reasonable statistical-mechanical 
model, for that matter — for any finite graph G. Rather, phase transitions arise 
only in the infinite-volume limit. That is, we consider some countably infinite graph 
Goo = (Voodoo) — usually a regular lattice, such as with nearest-neighbor edges 
- and an increasing sequence of finite subgraphs G n = (V n ,E n ). It can then be 

10 Indeed, some of the proofs become simpler when the multivariate generalizations are considered: 
see [74] for details. In particular, the multivariate approach clarifies the meaning of the mysterious 
number q = 32/27 [52,69]. 
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shown (under modest hypotheses on the G n ) that the (limiting) free energy per unit 
volume 

fGoo(q,v) = lim \Vn\~ 1 log Z Gn (q,v) (5.1) 

n—xx) 

exists for all nondegenerate physical values of the parameters 11 , namely either 

(a) q integer > 1 and — 1 < v < oo (using (2.7), see [68, Section 1.2]) 

or (b) q real > 0 and 0 < v < oo (using (1.1), see [61, Theorem 4.1] 

and [60,116]) 

This limit f Goo (q. v ) is in general a continuous function of v [and of q in case (b)]; but 
it can fail to be a real-analytic function of v and/or q. because complex singularities 
of log Z Gn (q , v) — namely, complex zeros of Z Gn (q , v) — can approach the real axis 
in the limit n —> oo. Therefore, the possible points of physical phase transitions 
are precisely the real limit points of such complex zeros. As a result, theorems that 
constrain the possible location of complex zeros of the partition function are of great 
interest. 12 

In particular, if it can be proven that the partition function Z Gn is nonvanishing 
in some complex domain D that does not vary with the size of the graph G n , then 
(under mild hypotheses) the infinite-volume free energy f Goo will be analytic in D 
— that is, in physical terms, there will be no phase transitions in D. Study of the 
complex zeros of the partition function is thus one way of proving the absence of 
phase transitions in specified regions of parameter space. 

It is useful to distinguish between “soft” and “hard” theorems of this type. 
The “soft” theorems assert the nonvanishing of the partition function Zq for some 
specified statistical-mechanical model in some “natural” large domain of multidi¬ 
mensional complex space (e.g. the unit polydisc or the product of right half-planes) 
for arbitrary graphs G. The “hard” theorems, by contrast, assert the nonvanishing 
of Zq in a domain that depends quantitatively on some parameter associated to G , 
such as its maximum degree. (“Hard” theorems are of physical interest only when 
the parameter in question is “local”, so that the domain obtained is uniform in the 
G n . E.g. maximum degree is OK, but total number of vertices is not.) 

The “soft” theorems are particularly interesting from a mathematical point of 
view: they assert that some combinatorially natural generating polynomial is non- 

11 Here “physical” means that the weights are nonnegative, so that the model has a probabilistic 
interpretation; and “nondegenerate” means that we exclude the limiting cases v = — 1 in (a) and 
q = 0 in (b), which cause difficulties due to the existence of configurations having zero weight. 

12 I don’t want to give the impression that the partition function is the only quantity of interest 
to statistical physicists. Indeed, what is really of physical significance is not the partition function, 
but rather the Boltzmann-Gibbs probability measure, which gives the probabilities for different 
configurations of the system: in the Potts (resp. random-cluster) model this measure is given by 
the summand of (2.7) [resp. (2.10)] divided by the partition function Zdq, v). From this point 
of view, the partition function is merely a normalization factor, of no intrinsic physical interest 
whatsoever! Nevertheless, derivatives of the partition function with respect to the parameters 
occurring in it yield correlation functions (i.e., expectations of particular random variables with 
respect to the Boltzmann-Gibbs measure), which are of direct physical interest. Furthermore, the 
complex zeros of the partition function are of interest for the reason just explained; the Yang-Lee 
picture thus provides one method for investigating phase transitions, to be used in conjunction with 
other approaches. 
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vanishing in some large domain of C n . I am aware of only three classes of such 
theorems (but I urge readers to try to discover new ones!): 

1) The Lee-Yang theorem for the ferromagnetic Ising model at complex magnetic 
field [83] and its generalizations. 13 

2) The Heilmann-Lieb theorem for the matching polynomial [64], 14 

3) The half-plane theorem for the spanning-tree polynomial and its generaliza¬ 
tions [36]. 

The Lee-Yang and Heilmann-Lieb theorems can, in fact, be given a unified combina¬ 
torial formulation: Let us equip the graph G = ( V. , E) with edge weights A = {A e } e g_E 
and vertex weights t = {L,}ievb and let us form the generating polynomial 

p G (t,A) = Wiim (n-W) m 

ACE KeCA J \i£V J 

where 


Lee-Yang: 


Wi{A) 


1 if deg^(i) is even 
ti if deg^(f) is odd 


(5.3) 


1 if deg A (i) = 0 

Heilmann-Lieb: Wi(A ) = t t if deg j4 (i) = 1 (5-4) 

.0 if deg A (i) > 1 

Then the Lee-Yang and Heilmann-Lieb theorems assert that Pq{ t, A) ^ 0 whenever 
A e > 0 for all e and Re t, t > 0 for all i (or R eti < 0 for all i). Indeed, the Lee-Yang 
and Heilmann-Lieb theorems can be given a unified proof. 15 The half-plane theorem 
for the spanning-tree polynomial will be discussed in Section 6 below. 

For a prototypical example of a “hard” theorem, consider the multivariate gen¬ 
erating polynomial for independent sets of vertices in a graph G = (V, E): ie 

z g( w) = yi n Wi ’ ( 5 - 5 ) 

ACV *eA 
A independent 

13 A partial bibliography (through 1980) of generalizations of the Lee-Yang theorem can be found 
in [85]. 

14 This theorem is most often quoted in the form of its univariate corollary, namely that the roots 
of the univariate matching polynomial are all (negative) real. But it is the multivariate theorem 
that is truly fundamental. Two simple proofs of the multivariate Heilmann-Lieb theorem can be 
found in [36, Theorem 10.1 and Remark 4 following it]. 

15 See [36, Sections 4.7 and 4.8, and Remark 4 in Section 10.1]. See also [110, 111] for related 
results, proven by a different method. 

16 We recall that a set ACV is called independent (or stable ) if it does not contain any pair 
of adjacent vertices. This is totally unrelated to the concept of independent sets of elements in a 
matroid; it is unfortunate that the same word is used. 
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where w = {wi}i£v are (complex) vertex weights. [In statistical mechanics, (5.5) is 
called the grand partition function of the hard-core lattice gas on G; the weights w t 
are called activities or fugacities .] It can then be shown [114] that if G is a graph 
of maximum degree A, and \w%\ < (A — 1) a_1 /A a for all i E V, then Zq( w) ^ 0. 
Indeed, this is a corollary of a more general result [114] that provides a sufficient 
condition on a set of radii R = for Zq to be nonvanishing in the closed 

poly disc \wi\ < Ri . 1T 

These theorems belong, in fact, to a long line of results in the mathematical- 
physics literature asserting that the grand partition function of a repulsive lattice 
gas is nonvanishing in some polydisc at small fugacity. (Physically this corresponds 
to the absence of phase transition for a gas at low density.) Two main methods of 
proof are used: 

1) The coefficients in the Taylor expansion of log Zq( w) are interpreted combina- 
torially: this is the so-called Mayer expansion [132]. Then, by direct estimation 
involving some nontrivial combinatorics, the Mayer expansion is proven to be 
convergent in a suitable polydisc |iuj| < Ri [12,21,24-26,29,80,91,101,103, 
104, 115, 120]. It immediately follows that Zq is nonvanishing in this same 
polydisc. 

2) The nonvanishing of Zq in the closed polydisc \w t \ < Rj is proven directly, by 
induction on the number of vertices in G [45,46,114,121]. It then immediately 
follows that the Taylor series for log Zq is convergent in the interior of this 
polydisc. 

The foregoing results have a wider application than one might at first think, 
because the hard-core lattice gas is not merely one interesting statistical-mechanical 
model; it is, in fact, the universal statistical-mechanical model in the sense that any 
statistical-mechanical model living on a vertex set Vo can be mapped onto a gas 
of nonoverlapping “polymers” on Vo, i.e. a hard-core lattice gas on the intersection 
graph of Vo [120, Section 5.7]. 18 This construction, which is termed the “polymer 
expansion” or “cluster expansion”, is an important tool in mathematical statistical 
mechanics [11,21,25,59,115]; it is widely employed to prove the absence of phase 
transition at high temperature, low temperature, large magnetic held, low density, 
or weak nonlinear coupling. In the antiferromagnetic Potts model it can be used to 
prove the nonvanishing of Zc(q, v) in a region |g| > R (see Section 9). 

6 q = 0: Spanning trees, electric circuits, and all that 

As explained in Section 2.3, the spanning-tree polynomial Tq(v) can be obtained 
from the multivariate Tutte polynomial Zc(q, v) via the double limit in which q —> 0 
and v is infinitesimal. In this section we shall study some remarkable properties of 
the spanning-tree polynomial and its generalization to matroids. 

17 These results are, in turn, closely related to the Lovasz local lemma [55,56,124,125] in proba¬ 
bilistic combinatorics; see [114] for a detailed development of this connection. 

18 The intersection graph of a finite set S is the graph whose vertices are the nonempty subsets 
of S, and whose edges are the pairs with nonempty intersection. 
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6.1 The matrix-tree theorem 

Let G = ( "V,E) be a connected graph, and let L G (x) be the Laplacian matrix 
for G with edge weights x = {x e } e ^E'- 

L G (x) = (6.1) 

{Ek&22e~ik X e 

where e ~ ij denotes that e connects i to j. [Equivalently, L G (x) = BXB T where 
B is the directed vertex-edge incidence matrix for any orientation of G, and X = 
diag(x).] By construction, the row (and column) sums of L G (x) are all zero, so its 
determinant is zero. Now let i be any vertex of G, and let L G (x)\j be the matrix 
obtained from L G (x) by deleting the ith row and column. Kirchhoff [77] proved in 
1847 the following striking result (see also [15,34,35,92,94,146]): 


Proposition 6.1 (matrix-tree theorem) det L G (x)w is independent of i and 
equals T G (x), the generating polynomial of spanning trees in G. 

Many different proofs of the matrix-tree theorem are now available; one simple proof 
is based on the Cauchy- Binet theorem in matrix theory (see e.g. [92]). 

More generally, it turns out that each minor of L G (x) enumerates a suitable 
class of rooted spanning forests [1,33,93]. To formulate this result, suppose that 
I,J CP, and let us denote by Lg( x )\/.j the matrix obtained from L G (x) by deleting 
the columns / and the rows J; when I = J, we write simply Lq(x)\j. Then the 
“principal-minors matrix-tree theorem” states that 

detL G (x)\{ ilv .. jir } = ^2 \\ x e, (6-2) 

eE.F 

where the sum runs over all spanning forests F in G composed of r disjoint trees, 
each of which contains exactly one of the “root” vertices ii,... ,i r . This theorem 
can easily be derived by applying Proposition 6.1 to the graph in which the vertices 
ii,... ,i r are contracted to a single vertex. Furthermore, the “all-minors matrix-tree 
theorem” (whose proof is a bit more intricate, see [1, 33, 93]) states that for any 
subsets I, J of the same cardinality r, we have 

detL G (x)\ J>J = ^ t(F,I,J) Y[x e , (6.3) 

FeT(I\J) eeF 

where the sum runs over all spanning forests F in G composed of r disjoint trees, 
each of which contains exactly one vertex from I and exactly one vertex (possibly 
the same one) from J; here e(F, /, J) = ±1 are signs whose precise definition is not 
needed here. 

The virtue of the matrix-tree theorem is that enumerative questions about span¬ 
ning trees (and rooted spanning forests) can be reduced to linear algebra. 
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6.2 Electric circuits and the half-plane property for graphs 

Now let us consider the graph G as an electrical network: to each edge e we 
associate a complex number x e , called its conductance (or admittance). 19 [The 
resistance (or impedance) is l/x e .] Suppose that we inject currents J = {Ji}iev 
into the vertices. What node voltages ip = {ifi}i & v will be produced? Applying 
Kirchhoff’s law of current conservation at each vertex and Ohm’s law on each edge, 
it is not hard to see that the node voltages and current inflows satisfy the linear 
system 

L g (x)<£> = J . (6.4) 

It is then natural to ask: Under what conditions does this system have a (unique) 
solution? Two obvious constraints arise from the fact that the row and column sums 
of Lq(x) are zero: firstly, the current vector must satisfy Yliev = 0 (“conservation 
of total current”), or else no solution will exist; and secondly, if <p is any solution, 
then so is <p + cl for any c ( “only voltage differences are physically observable” ) . 
So let us assume that Yliev = an d ^ us break the redundancy in the solution 
by fixing the voltage to be zero at some chosen reference node iq £ V (“ground”). 
Does the modified system 

-^G( x )\i 0 J\io (b - b) 

then have a unique solution? This will be so if and only if det Lq(x)\ iq is nonzero — 
which, by virtue of the matrix-tree theorem, is equivalent to T G (x) being nonzero. 

Simple counterexamples show that this is not always the case. Suppose, for 
instance, that G consists of a pair of vertices connected by two edges e, / in parallel. 
With x e = 1 and Xf = —1, it is easy to see that no solution exists (except when 
J = 0). Of course, the reader will object that negative resistances are physically 
unrealizable! Fine: consider instead x e = i and Xf = —i; once again there is 
no solution (unless J = 0). This example corresponds physically to a perfectly 
lossless capacitor together with a perfectly lossless inductor, exactly at their resonant 
frequency. But once again, the reader will rightly object that perfectly lossless 
components are physically unrealizable; every component in the real world exhibits 
some dissipation. This reasoning leads us to conjecture, on physical grounds, that 
if Rex e > 0 for all e (i.e. every branch is strictly dissipative), then the network is 
uniquely solvable once we fix the voltage at a single reference node iq € V. This 
conjecture turns out to be true, and we have: 

Theorem 6.2 Let G be a connected graph. Then the spanning-tree polynomial Tq 
has the “half-plane property” (HPP), i.e. Rex e > 0 for all e implies Tq(x) f 1 0. 

The proof of Theorem 6.2 is not difficult: Consider any nonzero complex vector 
ip = satisfying p lQ = 0. Because G is connected, we have B T p / 0. 

Therefore, the quantity 

P*L g (x)p = p*BXB T p = ^|(R T ^) e | 2 x e (6.6) 

e€E 

19 Complex admittances are relevant, for instance, in alternating-current circuits. Indeed, the 
reader is invited to imagine that we are considering an alternating-current circuit at some fixed 
frequency w. 
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has strictly positive real part whenever Re x e > 0 for all e; so in particular (BXB T ip)i ^ 
0 for some i ^ io- It follows that the submatrix of Tg( x ) obtained by suppressing the 
?'oth row and column is nonsingular, and so has a nonzero determinant. Theorem 6.2 
then follows from the matrix-tree theorem. 20 

It cannot be overemphasized how remarkable this theorem is. Suppose, for in¬ 
stance, that G is the complete graph on n vertices; then Tg( x ) is a homogeneous 
polynomial of degree n — 1 in n(n — l)/2 variables, containing n n ~ 2 monomials. It 
seems utterly miraculous, at first sight, that a polynomial of such complexity should 
be nonvanishing whenever Re x e > 0 for all e. The fact that it is so clearly expresses 
a deep property arising from the underlying combinatorial structure. 

An immediate corollary of Theorem 6.2 is that the complementary spanning-tree 
polynomial 



also has the half-plane property, since the map x e i—> l/x e takes the right half-plane 
onto itself. 

6.3 The half-plane property for matroids 

From a combinatorial point of view, the noteworthy fact is that the spanning 
trees of G constitute the bases of the graphic matroid M(G ), and their comple¬ 
ments constitute the bases of the cographic matroid M*(G). So Tq and Tq are the 
multivariate basis generating polynomials for M(G ) and M*(G), respectively. This 
naturally suggests generalizing Theorem 6.2 to more general matroids and, perhaps, 
to more general set systems. Before posing these questions precisely, we need to fix 
some notation and terminology. 

A set system (or hypergraph ) S on the (finite) ground set E is simply a collection 
S of subsets of E. Given any set system S on E, we define its (multivariate) 
generating polynomial to be 

p s{x) = Tl Xe ( 6 - 8 ) 

S£S e&S 

The rank of a set system is the maximum cardinality of its members (by convention 
we set rank = — oo if 5 = 0); equivalently, it is the degree of the generating poly¬ 
nomial Pg. A set system S is r-uniform if \S\ = r for all S € S, or equivalently if 
its generating polynomial Pg is homogeneous of degree r. We shall be particularly 
interested in the basis generating polynomial of a matroid M , 

p B(M)(x) = B M (x) = ^2 \\ x e ( 6 - 9 ) 

BeB(M) e£B 

We can now pose the following questions concerning possible extensions of The¬ 
orem 6.2: 

20 This proof is well known in the circuit-theory literature: see e.g. [35, Section 2.7] as well as 
the related results in [42, pp. 398-401, 430-431 and 850-851] [100, pp. 52-53 and 67-69]. It has, 
moreover, a natural physical interpretation: if g> = are the node voltages, then the real 

part of the quadratic form (6.6) is the total power dissipated in the circuit. 
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Question 6.3 For which matroids M does the basis generating polynomial Pbim) 
have the half-plane property? 

More generally: 

Question 6.4 For which r-uniform set systems S does the generating polynomial 
Ps have the half-plane property? 

These questions were recently studied in a long paper by Choe, Oxley, Wagner 
and myself [36]. Our original conjecture was that all matroids (and no non-matroidal 
set systems) have the half-plane property. That would be nice and neat, but it 
turns out to be false; and the truth is considerably more interesting and subtle. 
Our conjecture is half right: an r-uniform set system with the half-plane property is 
necessarily the set of bases of a matroid [36, Theorem 7.1]. But not every matroid 
has the half-plane property, and we do not yet have a complete characterization of 
those that do. Nevertheless, we can find large classes of matroids with the half-plane 
property: 

(a) Every complex unimodular matroid — or what is equivalent, sixth-root-of- 
unity matroid [137] — has the half-plane property [36, Theorem 8.1 and 
Corollary 8.2], 21 This class properly includes the regular matroids, which 
in turn properly includes the graphic and cographic matroids. The proof of 
the half-plane property for complex unimodular matroids is, in fact, a direct 
generalization of the proof just given for Theorem 6.2. 

(b) Every uniform matroid has the half-plane property. (Indeed, every loopless 
uniform matroid has the Brown-Colbourn property, which is stronger than 
the half-plane property: see Section 7 below.) 

(c) A significant subclass of transversal matroids — those we call “nice” — have 
the half-plane property [36, Section 10]. 22 

(d) All matroids of rank or corank at most 2 have the half-plane property [36, 
Corollary 5.5], as do all matroids on a ground set of at most 6 elements [36, 
Proposition 10.4], 

(e) The class of matroids with the half-plane property is closed under minors, 
duality, direct sums, 2-sums, series and parallel connection, full-rank matroid 
union, and some special cases of principal truncation, principal extension, prin¬ 
cipal cotruncation and principal coextension [36, Section 4]. 

Furthermore, we show that: 

21 It is proven in [36, Theorem 8.9] that a matroid is complex unimodular if and only if it is 
sixth-root-of-unity. 

22 In [36, Conjecture 13.16 and Question 13.17] it was conjectured that all rank-3 transversal 
matroids have the half-plane property, and the question was raised whether it might even be true 
that all transversal matroids have the half-plane property. Both conjectures are in fact false. Choe 
and Wagner [37] have exhibited a rank-4 transversal matroid that lacks the half-plane property (and 
indeed lacks the weaker Rayleigh property, see Section 6.4 below). I have found a class of rank-3 
transversal matroids that lack the half-plane property, and plan to publish them elsewhere [123]. 
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(f) A binary matroid has the half-plane property if and only if it is regular [36, 
Corollary 8.16]. 

(g) A ternary matroid has the half-plane property if and only if it is a sixth-root- 
of-unity matroid [36, Corollary 8.17]. 

Finally, we can show that certain matroids do not have the half-plane property: 
among these are the Fano matroid F 7 . the non-Fano matroid Ff , their relaxations 
F 7 , F 7 -3 and M(K±) + e, the matroids Pg, Pg and P", the Pappus and non-Pappus 
matroids, the free extension (non-Pappus \ 9) + e, and all their duals [36, Section 11], 
The first six of these examples are minor-minimal, and we conjecture that the others 
are as well; but we strongly suspect that this list is incomplete, and indeed we 
consider it likely that the set of minor-minimal non-half-plane-property matroids is 
infinite. 

More generally, we consider homogeneous multiaffine polynomials 

P( x ) = ^2 a s n^e (6.10) 

SCE,\S\=r ees 

with arbitrary complex coefficients as (not necessarily 0 or 1). We prove two neces¬ 
sary conditions for P ^ 0 to have the half-plane property: 

(a) P must have the “same-phase property”, i.e. all the nonzero coefficients as 
must have the same phase [36, Theorem 6.1]. So without loss of generality we 
can assume that all the as are real and nonnegative. 

(b) The support supp(P) = {S C E: as 0} must be the collection of bases of a 
matroid [36, Theorem 7.1]. 

This latter fact is particularly striking: it shows that matroids arise naturally from a 
consideration of homogeneous multiaffine polynomials with the half-plane property. 
We do not know whether the converse of this result is true, i.e. whether for every 
matroid M there exists a homogeneous multiaffine polynomial P with the half¬ 
plane property such that supp(P) = B{M). But it is true, at least, for all matroids 
representable over C [36, Corollary 8.2], 

We also also give two sufficient conditions for a homogeneous multiaffine poly¬ 
nomial P to have the half-plane property (or be identically zero): 

(a) Determinant condition [36, Theorem 8.1]: as = |det(A [ S)\ 2 for some r x n 
complex matrix A [here n = \E\, and A [ S denotes the square sub matrix 
of A using the columns indexed by the set S]. This corresponds to P(x) = 
det(AAA*) where X = diag(x) and * denotes Hermitian conjugate. 

(b) Permanent condition [36, Theorem 10.2]: as = per(A f S) for some r x n 
nonnegative matrix A. This corresponds to P(x) = per(AA). 

Unfortunately, the relationship between these sufficient conditions and the half¬ 
plane property looks complicated. Neither family of polynomials contains the other; 
their intersection is nonempty; and their union is a proper subset of the set of all 
homogeneous multiaffine polynomials with the half-plane property. 

In any case, the matroids with the half-plane property form a very natural class, 
which deserves further study. A long list of open questions can be found in [36, 
Section 13]. 
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6.4 The Rayleigh property for graphs and matroids 

Consider once again the graph G as an electrical network, with conductances 
{x e } e& E on the edges, and select a pair of distinct vertices i.j G V. Then Kirchhoff 
[77] showed in 1847 that the effective conductance between nodes i and j is 


Xj(x) 


Jb(x) 

T G /ij{x) ’ 


( 6 . 11 ) 


where G/ij is the graph obtained from G by contracting i and j to a single vertex. 23 

Suppose now (in contrast to the preceding subsections) that all the conductances 
x e are real and positive. Physical intuition suggests that increasing one or more of 
the branch conductances x e cannot cause the effective conductance between any pair 
of nodes to decrease ; in other words, we conjecture that Tjj(x) is a nondecreasing 
function of all the x e on [0, oo) E . Let us rephrase this statement by considering 
the graph H obtained from G by adjoining a new edge / connecting i to j; then 
G = H \ f and G/ij = H/f. Using 


T H\f(x) = Jfl\{ e j } (x) + x e T( H \ f y e (x) (6.12a) 

r H//(x) = T (/V/)\e(x) + x e T H/ {e ,/}(x) (6.12b) 


we hnd 


w e y -’ (x) 


T (H/e)\f (x) T [H/ f)\ e (x) - T H/{eJ] (x) T H \ {eJ} (x) 
T H/f(x ) 2 

^/e(x)T g //(x) - r g/{eJ} (x)T g (x) 

T H//(x) 2 


(6.13a) 

(6.13b) 


The conjecture then becomes: 


Theorem 6.5 (Rayleigh property for graphs) Let H = (V. E) be a finite undi¬ 
rected graph, and let e, f G E with e//. Then 

T H /e(x)T H /f(x) - T H /{ e jy (x) T h (x) > 0 (6.14) 


for all x > 0. 

Please note that (6.14) has a beautiful probabilistic interpretation: Consider the 
probability measure on spanning trees T C H giving weight (EU T*e)/T H (x) to 
the tree T. Then (6.14) asserts that the events e G T and / G T are negatively 
correlated. 


23 Proof: By definition, l/34,( x ) is the voltage difference induced between nodes i and j if we 
inject 1 ampere of current into i and extract 1 ampere from j. So let node j be “ground”, and 
apply (6.5) with io = j and J = <5; — 8j. We then have 


1 /Tij(x) = [LcMydii 1 


det lG(x)\{jj} 
det 1g(x)\j 


by Cramer’s rule. The principal-minors matrix-tree theorem (6.2) then tells us that det I/g(x)\j = 
T’g(x) and det L G (x\ {iij y = T a/ij (x). 
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Theorem 6.5 has both algebraic and variational/probabilistic proofs: see e.g. 
[4, Section 3.8] [57, Theorem 2.1] [37, Theorems 5.2 and 5.6] for the former and 
[48,87,88,99] for the latter. 

In an important recent paper, Choe and Wagner [37] have investigated the gen¬ 
eralization of this problem to matroids. Let M be a matroid with ground set E, 
and let Bm(x) = ^s(M)( x ) be its (multivariate) basis generating polynomial. We 
say that M is a Rayleigh matroid (or has the Rayleigh property ) in case 

%(x)%(x) - %/{e,/}(x) -Bm(x) > 0 (6.15) 

for all e, / € E with e/ / and all x > 0. Choe and Wagner [37] have shown that: 

(a) The class of Rayleigh matroids is closed under minors, duality, direct sums, 
and 2-sums. 

(b) A binary matroid is Rayleigh if and only if it does not contain S% as a minor. 
Equivalently, a binary matroid is Rayleigh if and only if it can be constructed 
from regular matroids, F-j. F j and AG(3, 2) by taking direct sums and 2-sums. 

(c) Every matroid with the half-plane property is Rayleigh (but the converse is 
false). 

(d) Every matroid on a ground set of at most 7 elements is Rayleigh. 

More recently, Wagner [134] has also shown that 

(e) Every matroid of rank or corank at most 3 is Rayleigh. 

On the other hand, Choe and Wagner [37] have exhibited a rank-4 transversal ma¬ 
troid that is not Rayleigh. 

The most striking of these results, in my opinion, are (b) and (c). For binary 
matroids, (b) provides a pair of beautiful characterizations of the Rayleigh property: 
one by excluded minors and the other by internal structure. It is at least conceivable 
that analogous results could be obtained for ternary matroids, GI ? (4)-representable 
matroids, or even arbitrary matroids. As for (c), one might have thought, a priori, 
that the half-plane property and the Rayleigh property are distinct but unrelated 
properties of electric circuits (one for complex edge conductances, the other for 
positive real edge conductances), which thus hold for all graphic matroids but extend 
to distinct and unrelated families of non-graphic matroids. Item (c) shows that the 
truth is quite different: the half-plane property is strictly stronger than the Rayleigh 
property. Indeed, in some vague sense it seems that the half-plane property is quite 
a bit stronger than the Rayleigh property: non-HPP matroids seem to be plentiful, 
while non-Rayleigh matroids seem to be fairly rare. In any case, the Rayleigh 
matroids form a very natural class and deserve further study; see [37] for a list 
of open problems. 

7 gr = 0 again: The reliability polynomial 
7.1 The Brown-Colbourn property for graphs 

Once again let G = {V, E) be a connected graph, considered now as a communi¬ 
cations network with unreliable communication channels: the edge e is assumed to 
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be operational with probability p e and failed with probability 1 — p e , independently 
for each edge. Let Rc(p) be the probability that every node is capable of commu¬ 
nicating with every other node (this is the so-called all-terminal reliability). Clearly 
we have 

Rg{ p) = W Pe II ( 1_ Pe), (7.1) 

ACE e£A e€E\A 
(V, A) connected 

where the sum runs over all connected spanning subgraphs of G , and we have written 
p = {p e } e£ E- We call Rq( p) the (multivariate) reliability polynomial [38] for the 
graph G ; it is a multiaffine polynomial, i.e. of degree at most 1 in each variable 
separately. If the edge probabilities p e are all set to the same value p, we write the 
corresponding univariate polynomial as Rg(p ), and call it the univariate reliability 
polynomial. We are interested in studying the zeros of these polynomials when the 
variables p e (or p) are taken to be complex numbers. 

The reliability polynomial Rg( p) is, of course, just the connected-spanning- 
subgraph polynomial C'g(v) in disguise: 

na-ft) 

.eCE 

Yl(l + y e) 

,e£E 

So we can equally well work with C'c(v), making the change of variables v e = 
Pe/(1 - Pe)- 

Brown and Colbourn [16] studied a number of examples and made the following 
conjecture: 

Univariate Brown—Colbourn conjecture. For any connected graph 
G, the zeros of the univariate reliability polynomial Rg(p ) all lie in the 
closed disc \p — 1| < 1 . In other words, if \p — 1| > 1 , then Rg{p) 7 ^ 0. 

Subsequently, I proposed [121] a multivariate extension of the Brown-Colbourn con¬ 
jecture: 

Multivariate Brown—Colbourn conjecture. For any connected graph 
G, if \p e — 1| > 1 for all edges e, then Rg(p) 0. 

In terms of Cq, the multivariate Brown-Colbourn conjecture states that if G is a 
loopless connected graph and |1 + v e \ < 1 for all edges e, then Cciy) 7^ 0. Loops 
must be excluded because a loop e multiplies Cg by a factor 1 + v e but leaves Rg 
unaffected. 

A few years ago, Wagner [133] proved, using an ingenious and complicated con¬ 
struction, that the univariate Brown-Colbourn conjecture holds for all series-parallel 
graphs. 24 As an illustration of the power of the “multivariate ideology”, let us now 

24 Unfortunately, there seems to be no completely standard definition of “series-parallel graph”; 
a plethora of slightly different definitions can be found in the literature [13,38,50,97,98]. So let 
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V 


1 + V 


(7.3) 



Rg( p) = 

Cg(v) = 
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prove the stronger result [121] that the multivariate Brown-Colbourn conjecture 
holds for all series-parallel graphs — and let us moreover do it in two lines: 2 ' 5 

Theorem 7.1 Let G be a loopless connected series-parallel graph. If |1 + v e \ < 1 
for all edges e, then Cq{w) / 0. 

Proof How should we prove a theorem for series-parallel graphs? The answer 
is obvious: just show that series and parallel reduction preserve the hypothesis 
|1 + v e \ < 1. If edges e\, e 2 (with weights v\, vf) are in parallel, they can be replaced 
by a single edge with weight u* satisfying 1 + u* = (1 + Ui)(l + U 2 ); obviously 
|l+Ui| < 1 for * = 1,2 implies that |l+u*| < 1. Likewise, if edges e\, e 2 (with weights 
^ 1 ,^ 2 ) are in series, they can be replaced by a single edge with weight u* satisfying 
1/u* = l/v\ + 1 /^ 2 , provided that we multiply Cq by the prefactor v\ + V 2 [cf. 
(4.23)/(4.24) specialized to q = 0]. Now |l + u| < 1 is equivalent to Re(l/u) < —1/2, 
so Re(l/uj) < —1/2 for i = 1,2 implies that Re(l/u*) < —1 < —1/2; furthermore, 
the prefactor v\ + V 2 in (4.24) is nonzero. Since the multivariate Brown-Colbourn 
conjecture manifestly holds for the base case of trees, it necessarily holds for all 
loopless connected series-parallel graphs. □ 

Does the multivariate Brown-Colbourn conjecture hold for all graphs? For sev¬ 
eral years I would have bet that it does (though I was unable to find a proof); but 
in 2002, Gordon Royle sent me the surprising news that the multivariate Brown- 
Colbourn conjecture fails already for the simplest non-series-parallel graph, namely 
the complete graph K 4 . The construction turns out to be very simple [108]. Since 
the univariate Brown-Colbourn conjecture holds for K 4 , let us try the next sim¬ 
plest situation, namely the bivariate one in which the six edges receive two different 
weights a and b. There are five cases: 

(a) One edge receives weight a and the other five receive weight b: 

C K4 (a,b ) = (86 3 + 56 4 + 6 5 ) + (8b 2 + 106 3 + 5b 4 + b 5 )a (7.4) 

(b) A pair of nonintersecting edges receive weight a and the other four edges receive 
weight b: 

C Ki (a,b ) = (46 3 + 6 4 ) + (86 2 + 86 3 + 26 4 )a + (46 + 66 2 + 46 3 + 6 4 )a 2 (7.5) 

(c) A pair of intersecting edges receive weight a and the other four edges receive 
weight b: 

CK 4 (a,b) = (36 s + 6 4 ) + (106 2 + 86 3 + 26 4 )a + (36 + 66 2 + 46 3 + 6 4 )a 2 (7.6) 


me be completely precise about my own usage: I shall call a loopless graph series-parallel if it 
can be obtained from a forest by a finite sequence of series and parallel extensions of edges (i.e. 
replacing an edge by two edges in series or two edges in parallel). I shall call a general graph 
(allowing loops) series-parallel if its underlying loopless graph is series-parallel. Some authors write 
“obtained from a tree”, “obtained from AV’ or “obtained from C 2 ” in place of “obtained from a 
forest”; in my terminology these definitions yield, respectively, all connected series-parallel graphs, 
all connected series-parallel graphs whose blocks form a path, or all 2-connected series-parallel 
graphs. See [13, Section 11.2] for a more extensive bibliography. 

25 This proof can be found in [121, Remark 3 in Section 4.1] or in [108, Theorem 5.6(c) =>■ (a)]. 
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(d) A 3-star receives weight a and the complementary triangle receives weight b: 

Ca' 4 (o, b) = (9b 2 + 36 3 )a + (66 + 9b 2 + 3b 3 )a 2 + (1 + 36 + 3b 2 + b 3 )a 3 (7.7) 

(e) A three-edge path receives weight a and the complementary three-edge path 
receives weight b: 

C K4 (a, b) = b 3 + (7b 2 + 3b 3 )a + (7b + 9b 2 + 36 3 )a 2 + (1 + 36 + 36 2 + b 3 )a 3 (7.8) 

Now plot the roots a when b traces out the circle |1 + b\ = 1, or vice versa. In cases 
(b) and (d) it turns out that the roots can enter the “forbidden discs” |1 + a\ < 1 
and |1 + b\ < 1 (see [108, Figures 1 and 2]). In fact, this is quite easy to understand 
analytically: if we solve the equation Ck 4 (ci, b) = 0 for a in terms of b , expanding in 
power series for b near 0, in cases (b) and (d) we obtain 

a = 5 l b + 5 2 b 3/2 + 0(b 2 ) (7.9) 

with di < 0 and 5 2 ^ 0. If we now set b = — 1 + e , then one of the roots will have 

I + a| < 1 for small 6^0 (arising from the S 2 b 3 / 2 term). A small perturbation (so 
that |1 + b\ < 1) then yields a counterexample to the multivariate Brown-Colbourn 
conjecture. 

In fact, with a little more work Royle and I are able to prove the following [108, 
Theorem 5.6]: 

Theorem 7.2 A loopless connected graph G has the multivariate Brown-Colbourn 
property if and only if it is series-parallel. 

Moreover, as a corollary of the failure of the multivariate Brown-Colbourn con¬ 
jecture for the complete graph K t i, Royle and I are able to show [108] that the 
univariate conjecture is false as well: counterexamples include a 4-vertex, 16-edge 
planar graph that can be obtained from K t ± by adding parallel edges, and a 1512- 
vertex, 3016-edge simple planar graph that can be obtained from K 4 by adding 
parallel edges and then subdividing edges. These counterexamples to the univariate 
conjecture are fairly easy to find once one has in hand the K 4 counterexample to 
the multivariate conjecture, but would probably not otherwise be guessed. This 
illustrates once again the advantages of considering the multivariate problem, even 
if one is ultimately interested in a particular univariate specialization. 

7.2 The Brown-Colbourn property for matroids 

The foregoing considerations can be extended to matroids. Let M be a matroid 
with ground set E, and let Sm (v) be its spanning-set polynomial. We then say that 
M has the (multivariate) Brown-Colbourn property in case Sm (v) / 0 whenever 

II + v e | < 1 for all e G E , 26 

Gordon Royle and I are currently investigating the Brown-Colbourn property 
for matroids [109]. Here are some of our results so far: 

26 Unfortunately, in [36] we used the dual definition (i.e., a matroid M has the Brown-Colbourn 
property in the sense of [36] if and only if M* has the Brown-Colbourn property as defined here). 
I think the present definition is more natural in view of the relation to graphs (namely, G has the 
Brown-Colbourn property as defined here if and only if its cycle matroid M(G) does). I therefore 
pledge to stick to this definition in the future. 
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(a) The class of matroids with the Brown-Colbourn property is closed under dele¬ 
tion, direct sums, 2-sums, series connection and parallel connection. 2 ' 

(b) If M has the Brown-Colbourn property, then it has the half-plane property as 
well [36]. But the converse does not hold, as the example of M{K t i) shows. 

(c) If M has the Brown-Colbourn property, it must be loopless. 

(d) Every loopless uniform matroid (i.e., every uniform matroid U r , n with 1 < r < 
n) has the Brown-Colbourn property [36,133]. 

(e) A binary matroid has the Brown-Colbourn property if and only if it is the 
cycle matroid of a loopless series-parallel graph. 

Most of these results are fairly easy to prove. 

8 q = 0 yet again: Spanning forests 

As discussed in Section 6.1, Kirchhoff’s matrix-tree theorem [15,34,77,92,94,146] 
and its generalizations to arbitrary minors [1,33,93] express the generating polyno¬ 
mials of spanning trees and rooted spanning forests in a graph G as determinants 
associated to the graph’s Laplacian matrix. Here I would like to discuss a recently 
discovered extension of the matrix-tree theorem that provides a compact represen¬ 
tation also for the generating polynomial of unrooted spanning forests, Fq{ w). 

Recall first that one useful formalism for manipulating determinants is Grass- 
mann algebra; indeed, any determinant can be written as a Gaussian “integral” over 
Grassmann variables. 28 In particular, it follows from the matrix-tree theorem that 
the generating polynomials of spanning trees and rooted spanning forests can be 
written as Gaussian Grassmann integrals. 

Very recently, Caracciolo, Jacobsen, Saleur, Sportiello and I [30] have proven 
a generalization of the matrix-tree theorem in which a large class of combinatorial 
objects are represented by suitable non-Gaussian Grassmann integrals. As a special 
case, we show that the generating polynomial of unrooted spanning forests, which 
arises as a q —> 0 limit of the multivariate Tutte polynomial [cf. (2.18)], can be 
represented by a Grassmann integral involving a Gaussian term and a particular 
quartic term. Although this representation has not yet led to any new rigorous 
results concerning w), it has led to important non-rigorous insights into the 
behavior of Fq( w) for large subgraphs of a regular two-dimensional lattice [30,75] 

insights that may yet be translatable into theorems by exploiting the rigorous 
renormalization-group methods developed in recent decades by mathematical physi¬ 
cists (see e.g. [107]). 

As in Section 6.1, let G = (V. E) be a finite undirected graph, and let L = 
Lg{ w) be the Laplacian matrix for G with edge weights w = {w e } e& E- Let us now 

27 The class of matroids with the B-C property fails to be closed under contraction, at least for a 
trivial reason: namely, contraction of one of a set of parallel elements will produce one or more loops, 
which are incompatible with the B-C property [see item (c) below]. But we do not know whether 
the contraction of a simple matroid with the B-C property can fail to have the B-C property. 

28 For introductions to Grassmann algebra and Grassmann-Berezin integration, see e.g. [1,147]. 
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introduce, at each vertex i € V, a pair of Grassmann variables Ai A- All of these 
variables are nilpotent of order 2 (i.e., iff = Al = 0), anticommute, and obey the 
usual rules for Grassmann integration [1,147]. Writing P(V’, V’) = YlieV^A^A, we 
have, for any matrix A, 

Jv ( ^)e^ = d etA (8.1) 

and more generally 

J$) A^jl • • • ArAr 

e(ii,..., |ji,..., jV) det ,...,ir-},{yi 

where the sign e(*i,..., i r \ji, ■ ■ • ,j r ) = ±1 depends on how the vertices are ordered 
but is always +1 when (?4 ,... ,i r ) = (j i,... ,jr )• These formulae allow us to rewrite 
the matrix-tree theorems in Grassmann form; for instance, (6.2) becomes 

/ w, vo (n V’^V’iJ = yi n We • ( 8 - 3 ) 

J \a=l / F£F(h,...,i r ) e£F 

Let us next introduce, for each connected (not necessarily spanning) subgraph 
T = (Vr,-Bp) °f G, the object 


Q r 



(8.4) 


(Note that each Q p is even and hence commutes with the entire Grassmann algebra.) 
Now consider an unordered family T = {Ti,... ,T/} with l > 0, and let us try to 
evaluate an expression of the form 

jV(A$)Qr 1 ---Qr l e* L * . (8.5) 

If the subgraphs Ti, ... , T/ have one or more vertices in common, then this integral 
vanishes on account of the nilpotency of the Grassmann variables. If, by contrast, the 
Ti, ..., T; are vertex-disjoint, then (8.3) expresses JV(ip, VO (ill-i n t ey r AA^j 

as a sum over forests rooted at the vertices of Vp = Ul=i Tr k • I n particular, all the 
edges of £p = (J fc=1 k must be absent from these forests, since otherwise two 
or more of the root vertices would lie in the same component (or one of the root 
vertices would be connected to itself by a loop edge). On the other hand, by adjoin¬ 
ing the edges of £p, these forests can be put into one-to-one correspondence with 
what we shall call T-forests, namely, spanning subgraphs H in G whose edge set 
contains £p and which, after deletion of the edges in £p, leaves a forest in which 
each tree component contains exactly one vertex from Vp- (Equivalently, a T-forest 
is a subgraph H with l connected components in which each component contains 
exactly one T *, and which does not contain any cycles other than those lying entirely 
within the T,;. Note, in particular, that a T-forest is a forest if and only if all the 



The multivariate Tutte polynomial 


37 


Tj are trees.) Furthermore, adjoining the edges of Ep provides precisely the factor 
WeeE T We - Therefore 


Jv(^) Qr t ■ • • Qr t e^ 


e 

H&T r e£H 


( 8 . 6 ) 


where the sum runs over all T-forests H. 

We can now combine all the formulae (8.6) into a single generating function, by 
introducing a variable tr for each connected subgraph T of G. Since 1+trQr = e fr< 3r ; 
we have 

r 4>Lip+J2 trQr x r \ x r 

V^,^) e r = Y (II tr ) EE n^- ( 8 - 7 ) 

r vertex- Ter HeT T eeH 

disjoint 

We can express this in another way by interchanging the summations over T and H. 
Consider an arbitrary spanning subgraph H with connected components H \,..., i7); 
let us say that T marks Hi (denoted T -< Hi) in case Hi contains T and contains no 
cycles other than those lying entirely within T. Define the weight 

W{Hi) = Y *r • (8.8) 

r <Hi 

Then saying that H is a T-forest is equivalent to saying that each of its components 
is marked by exactly one of the Tj; summing over the possible families T, we obtain 

f _ V>Z/!/ , +]C trQr f J-r \ -r-r 

V^,^)e r = Y n^) m-e. (8.9) 

H spanning C G \*=1 / e£H 


This is our general combinatorial formula. Extensions allowing prefactors 
are also easily derived. 

Now consider the special case in which tr = t whenever T consists of a single 
vertex with no edges, tr = u whenever T consists of two vertices linked by a single 
edge, and tr = 0 otherwise. We have 




exp 


ipL-i/j +1 y v # i + u y 
i (tj) 

= e n ( 8 .io) 

F G T M=1 / e£F 

F=(F 1 ,...,F l ) 


where the sum runs over spanning forests F in G with components F \,..., Ff, here 
Vp z and | Ef { \ are, respectively, the numbers of vertices and edges in the tree 
Fi. [The four-fermion term ^ can equivalently be written, us¬ 
ing nilpotency of the Grassmann variables, as —(uj 2) If u = 0, 

this formula represents vertex-weighted spanning forests as a determinant (“massive 
fermionic free field”) [6,51]. More interestingly, since {Vp^ — | Ep t \ = 1 for each tree 
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F t , we can take u = —t and obtain the generating function of unrooted spanning 
forests with a weight t for each component. This is furthermore equivalent to giving 
each edge e a weight w e /t, and then multiplying by an overall prefactor We 
have therefore proven: 


Proposition 8.1 Let G = (V, E) be a finite undirected graph, let L be the Lapla- 
cian matrix for G with edge weights w = {w e } e& E, and let Fq be the generating 
polynomial of spanning forests in G. Then 


T > ('ip,'ip) exp fiLfi + t^ifi'fi + ^ y; ifjifiLjjipjipj = t^^Fafw/t) . (8.11) 


hi 


This representation of unrooted spanning forests is the translation to generating 
functions and Grassmann variables of a little-known but important paper by Liu 
and Chow [86]. 

It is an open question whether there exists some analogue of Proposition 8.1 for 
the object dual to Fg(w), namely the generating polynomial Cg(v) of connected 
spanning subgraphs, when the graph G is non-planar. 

It would also be interesting to know whether the foregoing identities are in any 
way related to the forest-root formula of Brydges and Imbrie [22,23,65,66]. 


9 Absence of zeros at large |q| 

Combinatorialists have long been interested in the real or complex zeros of the 
chromatic polynomial Pc(q). A fair number of interesting theorems have by now 
been proven, notably concerning the real zeros, but even here a vast number of open 
problems remain; moreover, very little is known rigorously concerning the complex 
zeros. (See [70] for an excellent recent survey treating both real and complex zeros.) 

In this section I would like to discuss one aspect of the complex-zero problem 
where some progress has recently been made by exploiting methods from mathe¬ 
matical statistical mechanics: namely, the absence of zeros at large enough |g|. 

9.1 Bounds in terms of maximum degree and its relatives 

The bounds I want to discuss [121] apply, in fact, not only to the chromatic poly¬ 
nomial but to the multivariate Tutte polynomial Zq((j, v) throughout the “complex 
antiferromagnetic regime’' |1 + v e \ < 1. These bounds come in several variants, but 
here is a typical one: 

Theorem 9.1 Let G = (V, E) be a loopless finite undirected graph of maximum 
degree A. Suppose that the edge weights v = {v e } e ^E satisfy |l+n e | < 1 for all e. Let 
Wiax = max|v e |. Then all the zeros of Zciq,^) lie in the disc |q| < 7.963907u max A. 

e£E 

For the chromatic polynomial (v e = — 1 for all e) one obtains the immediate corollary: 

Corollary 9.2 Let G = (V, E) be a loopless finite undirected graph of maximum 
degree A. Then all the zeros of the chromatic polynomial Pc(q ) He in the disc 
|?| < 7.963907A. 
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Corollary 9.2 answers in the affirmative a question posed by Brenti, Royle and 
Wagner [14, Question 6.1], generalizing an earlier conjecture of Biggs, Damerell and 
Sands [7] limited to regular graphs. 

Of course, the constant 7.963907 is an artifact of the proof; it is presumably 
far from sharp. The linear dependence on A is, however, best possible, since the 
complete graph Ka+i has chromatic roots 0,1, 2,... , A. Furthermore (and perhaps 
surprisingly), the complete graph Ka+i is not the extremal graph for this prob¬ 
lem, and a bound |qj < A is not possible. In fact, a non-rigorous (but probably 
rigorizable) asymptotic analysis, confirmed by numerical calculations, shows [112] 
that the complete bipartite graph R'a.a has a chromatic root aA + o(A), where 
a = — 2/W{— 2/e) « 0.678345 + 1.447937i; here W denotes the principal branch of 
the Lambert W function (the inverse function of w *—> we w ) [39]. So the constant in 
Corollary 9.2 cannot be better than |a| ~ 1.598960. 

A complete proof of Theorem 9.1 and its variants can be found in [121] (see 
also [11]); here I would like simply to sketch the method, which I hope will be of 
wider use. 

It is immediately apparent that Theorem 9.1 and Corollary 9.2 are “hard” results 
in the sense of Section 5: the zero-free region depends on the maximum degree A. 
The proof of these results follows the plan sketched there: namely, the multivariate 
Tutte polynomial on the graph G = (V, E) is first mapped onto a gas of nonover¬ 
lapping “polymers” on V, i.e. a hard-core lattice gas on the intersection graph of V ; 
then this hard-core lattice gas is controlled by means of the theorems mentioned at 
the end of Section 5. 

The mapping to a polymer gas is really quite simple; indeed, it is nearly trivial. 
The definition (1.2) writes the multivariate Tutte polynomial Z G (q, v) as a sum over 
spanning subgraphs ( [V ., A). Let us perform this sum in two stages: First we sum over 
partitions of V into disjoint nonempty subsets Vi,..., 14 (k > 1 arbitrary); these 
subsets will correspond to the vertex sets of the connected components of (V, A). 
Then we sum over the ways of choosing edges within each component so as to make 
it connected. From (1.2) we see that a component V) = S will get a weight 

w{S) = q 1-151 ^2 W Ve (9-la) 

BCE eCB 
(S, B) connected 

= q 1 ^ 1 C G \ S (y) (9.1b) 

where G\S denotes the induced subgraph. Note, in particular, that (provided G is 
loopless) any set S of cardinality 1 gets weight w(S) = 1. So we need not consider 
such sets explicitly: it suffices to define the “polymers” to be the sets V) of cardinality 
2 or more; we then know that any vertex not covered by a polymer must be covered 
by a set Vi of cardinality 1. The polymers are therefore an arbitrary family (possibly 
empty) of disjoint sets S) C V of cardinality 2 or more. We have thus written 
Z G (q, v) as an independent-set polynomial (= hard-core lattice gas) of the form 
(5.5) for a graph G = (V, E) whose vertices are the subsets of V of cardinality 2 or 
more, and whose edges are the pairs with nonempty intersection. The fugacities are 
given by (9.1). 

The usefulness of this representation for proving the absence of zeros at large 
|<?| comes from the fact that the fugacities w(S) are all suppressed by powers of q -1 
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(since [S'! > 2), hence are small for large |</|. Furthermore, this suppression operates 
more strongly for larger polymers. This raises the hope that, when \q\ is sufficiently 
large, the model lies in a “low-fugacity” regime where the methods mentioned in 
Section 5 can be brought to bear. This is in fact the case. 

Of course, I have sloughed over one crucial point: the sum over B in (9.1a). 
It is here, not surprisingly, that the maximum degree A enters (along with u max ). 
What one needs, it turns out, is a bound on the number of connected subgraphs 
H C G containing a specified vertex x and having a specified number of vertices 
m = |<S'|. What is easier to prove, however, is a bound on the number of connected 
subgraphs H C G containing a specified vertex x and having a specified number of 
edges [121, Section 4.2] [72], The relation between the two is far from obvious, but it 
turns out that in the complex antiferromagnetic region |l + u e | < 1 (and only there!), 
a beautiful inequality due to Oliver Penrose [101] allows connected subgraphs to be 
bounded by frees . 29 The upshot is that the sum over B in (9.1a) can be controlled, 
with the result that w(S) satisfies a uniform bound that is exponentially decaying 
in \S\, provided that \q\ is large enough compared to A. This suffices to ensure the 
convergence of the “polymer expansion”, and hence the nonvanishing of Zc(q,v). 

I believe it should be possible to extend this result beyond the complex antifer¬ 
romagnetic regime, by bounding directly the number of connected subgraphs having 
a specified number of vertices [73]. But the resulting bound on \q\ will no longer 
be linear in A and u max ; rather, it will behave roughly like v max . And this is not 
simply an artifact of the method of proof: in the (/-state Potts ferromagnet (u > 0, 
q > 0 ) on the the simple hypercubic lattice T, d with nearest-neighbor edges (hence 
A = 2d), there is a first-order phase transition at vfiq) = q l ^ d [1 + 0(l/q)\ for all 
sufficiently large q [10,79,81,82,90]. As discussed in Section 5, this means that for 
v near vt{q) there will be complex zeros of Zc(q,v) for large subgraphs G C T* d . 

By an ad hoc (and aesthetically unsatisfying) method, Theorem 9.1 and Corol¬ 
lary 9.2 can be strengthened so as to bound the roots, not in terms of the maximum 
degree A, but in terms of the second-largest degree A 2 [121, Section 6 ]. For simplic¬ 
ity let me state only the result for chromatic polynomials: 


Corollary 9.3 Let G = ( V ., E) he a loopless finite undirected graph of second-largest 
degree A 2 . Then all the zeros of the chromatic polynomial Pc(q) He in the union of 
the discs |(/| < 7.963907A2 and \q — 1 | < 7.963907A2. In particular, they all lie in 
the disc |g| < 7.963907A2 + 1 . 

It should be stressed that “second-largest” in Corollary 9.3 cannot be replaced by 
“third-largest”. Indeed, the generalized theta graphs ©C.p) [18,122] constitute a 
family of planar (in fact, series-parallel) graphs in which all but two vertices have 
degree 2 , yet their chromatic roots are dense in the region {9 6 C: \q — 1| > 1} 
[ 122]. 30 So, with one large-degree vertex, the chromatic roots remain bounded; but 
with two, all hell can break loose. 

29 See [121, Section 4.1] and [114, Section 2.2] for further discussion. 

30 The graph consists of a pair of endvertices connected by p internally disjoint paths each 

of length s. 
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9.2 Bounds in terms of maxmaxflow? 

Alas, all of the foregoing bounds have a severe defect: they are unnatural, because 
the multivariate Tutte polynomial Zc(q, v) factorizes over blocks [cf. (4.2)], while A 
and A 2 can grow arbitrarily large when blocks are glued together at a cut vertex. As 
a consequence, the bounds obtained from Theorem 9.1 and its variants are sometimes 
very far from sharp: for instance, a tree can have arbitrarily large A and A 2 , but 
its chromatic roots are only 0 and l! Of course, this defect can be cured by the 
ad hoc technique of applying Theorem 9.1 ff. to each block of G rather than directly 
to G. But the deeper significance of this remark is that the quantities A and A 2 
are unnatural for studying the multivariate Tutte polynomial because they are not 
matroidal. 

One would like, therefore, to find a graph invariant that is smaller than A and 
A 2 , that “trivializes over blocks” (and ideally generalizes to matroids), and that is 
strong enough to bound the roots of Pg(q) and Zc(q,v)- I have a candidate: it 
is the maxmaxflow [72,119,121], defined as follows: For x,y € V with x j- y, let 
A g(x, y ) be the max flow from x to y: 


^G(x,y) 

= max of edge-disjoint paths from x to y 
= min of edges separating x from y 

(9.2a) 

(9.2b) 

The maxmaxflow A (G) 

is then defined by 



A(G) = max A o(x,y) . 
x, y € V 
a ' + y 

(9.3) 


[Note the contrast with the edge-connectivity, which is the minimum of A o(x,y) 
over x y] Clearly A c(x,y) < min [dcix), dciy)}, so that 

A (G) < A 2 (G) . (9.4) 

Furthermore, it is easy to see that maxmaxflow “trivializes over blocks” in the sense 
that A(G) = maxi<i<6 A(Gj) where Gi,... ,G(, are the blocks of G. 

An apparently very different quantity can be defined via cocycle bases. For X, Y 
disjoint subsets of V, let E(X,Y) denote the set of edges in G between X and Y. 
A cocycle of G is a set E{X,X C ) where X C V and X c = V \ X. It is well-known 
that the cocycles of G form a vector space over GF{ 2) with respect to symmetric 
difference; this is called the cocycle space of G. Let A(G) be the minmax cardinality 
of the cocycles in a basis, i.e. 


A(G) = minmax [Cl, (9.5) 

B CgB 

where the min runs over all bases B of the cocycle space of G , and the max runs 
over all the cocycles C in the basis B. Since one special class of cocycle bases 
consists of taking the stars C(x) = E({x}, {x} c ) for all but one of the vertices in 
each component of G, we clearly have 


A(G) < A 2 (G) . 


(9.6) 
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The relationship, if any, between maxmaxflow and cocycle bases is perhaps not 
obvious at first sight. But Bill Jackson and I have proven [72] that 

A (G) = A (G) . (9.7) 

The two definitions thus give dual approaches to the same quantity. 

All this theory, extends in fact, to the more general situation of a finite undi¬ 
rected graph G equipped with nonnegative real edge weights w = {w e } e& E- (In the 
application to the multivariate Tutte polynomial, one might want to take w e = \v e \ 
or something similar.) 

The ultimate goal of this work is to prove (or disprove!) strengthenings of The¬ 
orem 9.1 ff. in which maximum degree or second-largest degree are replaced by 
maxmaxflow. For instance, for chromatic polynomials one would like to prove the 
following: 

Conjecture 9.4 There exist universal constants C( A) < oo such that all the chro¬ 
matic roots (real or complex) of all loopless graphs of maxmaxflow A lie in the disc 
|qj < C(A). Indeed, I conjecture that C(A) can be taken to be linear in A. 

It is natural to try to prove Conjecture 9.4 by modifying the arguments of [121] 
so as to decompose a spanning subgraph of G into its 2-connected (rather than 
merely connected ) components. Such an argument will require, at a minimum, a 
good bound on the number of 2-connected subgraphs H C G containing a specified 
edge e and having a specified number of vertices or edges. Bill Jackson and I have 
found some bounds of this type [72], But this is only a first step, and the other 
obstacles in generalizing this proof — such as controlling the interaction between 
the 2-connected components — have yet to be overcome. 

It is worth mentioning that A can be defined also for binary matroids, general¬ 
izing the definition for graphs. But we do not yet know whether it is strong enough 
to give a good (i.e., singly exponential) bound on the number of 2-connected dele¬ 
tion minors containing a specified element e and having a specified size. Nor do we 
know whether it is useful in bounding the zeros of the matroid chromatic polynomial 

Z M {q, -!)• 


9.3 Some further questions 

There is no compelling reason, a priori , to limit attention in the complex q- 
plane to discs centered at the origin; many other types of regions are of interest. For 
instance, the following conjecture might conceivably be true: 

Conjecture 9.5 All the chromatic roots (real or complex) of all loopless graphs of 
maxmaxflow A lie in the half-plane Req < A. 

By factoring Pg(q)/q and using the fact that its roots are either real or occur in 
complex-conjugate pairs, one easily deduces that the truth of Conjecture 9.5 would 
imply: 

Conjecture 9.6 If G has maxmaxflow A, then Pg{q)/q and all its derivatives are 
nonnegative at q = A. [The same then holds also for Pg(q) and all its derivatives.] 
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Since the truth of Conjecture 9.6 would imply that Pg(q) has all nonnegative Taylor 
coefficients at q = A, and since Pg{q) is not identically zero (when G is loopless), 
Conjecture 9.6 would imply: 


Conjecture 9.7 If G is a loopless graph of maxmaxflow A, then Pg{q) > 0 for 
q > A. In particular, Pg{q) has no roots in the real interval (A, oo). 


Even the weaker versions of these conjectures with A replaced by A 2 or A are 
open! So one has the 3x3 matrix of conjectures 


9.5 => 9.5' =^> 9.5" 

4 4 A 

9.6 =^> 9.6' =^> 9.6" 

4 4 A 

9.7 =^> 9.7' =► 9.7" 


in which the strongest might well be true and the weakest is still an open question. 31 
It might be useful to seek counterexamples, either by systematic calculation on small 
graphs up to ~ 20 vertices or, perhaps more fruitfully, by constructing suitable 
infinite families of graphs. It goes without saying that I have no idea how to prove 
any of these conjectures, or even any compelling reason to believe they are true 
(other than my inability thus far to find counterexamples). 

We can pose these questions most generally as follows: Let Q be a class of finite 
graphs, and let V be a subset of the complex plane. Then we can ask about the sets 

Si(G, V) = IJ IJ {qeC: Z G (q,v) = 0} (9.8) 

Geg vev 

S 2 (G,V) = U U Z G (g,v) = 0} (9.9) 

G&g v: Ve€V Ve 

Among the interesting cases are the chromatic polynomials V = {—1}, the antifer¬ 
romagnetic regime V = [—1,0], and the complex antiferromagnetic regime V = A = 
{v <G C: |1 + u| < 1}. 

The theorems discussed in this section show that the sets S-ziGr-, -4) and S 2 {G' r , -A) 
are bounded , where Q r (resp. Q' r ) is the set of all graphs whose maximum degree (resp. 
second-largest degree) is less than or equal to r. But boundedness is a rather weak 
statement about a set; one would like to learn more about its location in the complex 
plane. 

In the other direction, I have recently shown [122] that if Q is the family of all 
generalized theta graphs 0( S ’A> [18,122], then Si(f?,{—1}) is dense in the region 
{(/ E C: \q — 1| > 1}. Moreover, if Q' denotes the family of joins of @( S ’P) with the 
complete graph K 2 , then S\(QuG', {—1}) is dense in the entire complex plane. 32 The 

31 STOP PRESS!!! Gordon Royle (private communication) lias found a family of counterexam¬ 
ples to Conjectures 9.5 and 9.5'. His graphs are planar (but not series-parallel) and all vertices but 
one have degree 3 (hence A = A 2 = 3). Included in his family are a 47-vertex graph with chromatic 
roots q « 3.0129950712±0.8089628639 i, a 95-vertex graph with chromatic roots q « 3.0536525915± 
0.7547530551 i, a 191-vertex graph with chromatic roots q « 3.07174237056 ± 0.7105232675 i, and 
a 383-vertex graph with chromatic roots q « 3.0766232972 ± 0.6746120243 i. 

32 The join of G and H , denoted G + H, is the graph obtained from the disjoint union of G and 
H by adding one edge connecting each pair of vertices x € V(G), y £ V{H). 
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idea of the construction in [122] is to exploit the reduction formulae of Section 4.6: 
suitably chosen 2-rooted subgraphs (in this case single edges) are concatenated in 
order to create a larger 2-rooted subgraph that has some desired value of the effective 
coupling u e ff. Such a construction is obviously inapplicable to 3-connected graphs, 
in which 2-rooted subgraphs of more than a single edge cannot occur. It is thus 
reasonable to ask: What is the closure of the set of chromatic roots of 3-connected 
graphs? The answer, alas, is still the whole complex plane [71]. To see this, consider 
the graphs ©( s ’?d + K n for fixed n and varying s,p: they are (n + 2)-connected, and 
their chromatic roots (taken together) are dense in the region |g — (n + 1)| > 1. In 
particular, considering both n and n + 2, the chromatic roots are dense in the whole 
complex plane. So even arbitrarily high connectedness does not, by itself, stop the 
chromatic roots from being dense in the whole complex plane. 

On the other hand, the graphs Q^’Pl + K n are non-planar whenever n > 1 and 
p > 3. This suggests posing a more restricted question: 

Question 9.8 What is the closure of the set of chromatic roots for 3-connected 
planar graphs? 

Here the answer may well be much smaller than CorC\{|g — 1] < 1}. But it 
will not be a bounded set [71], since the bipyramids C n + K -2 are 4-connected plane 
triangulations, but their chromatic roots are unbounded [106,117,118]. 

9.4 A final remark 

The ultimate goal of this research is to prove theorems of the type “such-and-such 
graph property (or conjunction of properties) implies such-and-such bound on the 
chromatic roots”. But there is a virtually unlimited number of possible assertions 
of this type; the hard part is to figure out which ones are true! 

There is thus, in my opinion, much room for numerical experiment as a means 
for obtaining intuition about which graph properties (e.g. maximum degree and its 
relatives, girth, planarity, connectivity, ..) affect which aspects of the chromatic 
roots (e.g. large positive real roots, large positive or negative real part, large imag¬ 
inary part, etc.). Some information can be obtained by systematic calculation on 
all graphs of a specified type up to as many vertices as can be handled — or, for 
suitable classes (e.g. cubic graphs), by random generation of typical graphs with a 
rather larger number of vertices — in an attempt to correlate properties of the graph 
with properties of the chromatic roots. Guided by this information, one may then 
choose to focus on special classes of graphs that potentially exhibit some type of 
“extremal” behavior, and either pursue the numerical calculation to larger graphs 
or, better yet, find infinite families of graphs that exhibit such “extremal” behavior 
and for which one can compute exactly the asymptotic behavior of the chromatic 
roots. 

It is important to keep in mind that many asymptotic properties of chromatic 
roots are achieved very slowly as the number of vertices grows: for instance, the 
roots of cubic graphs move quite slowly towards negative real part [105], and the 
roots of the generalized theta graphs Q( s 'PI fill out the complex plane extremely 
slowly as s,p —> oo [18]. As a result, numerical experiments can be deceptive unless 
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interpreted with extreme care. It is for this reason that exact computations with 
well-chosen infinite families of graphs are particularly valuable. 
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